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D ■ 

C/} ■ Abstract. 

This paper deals with the introduction of a decomposition of the deformations of curved thin beams, 
with section of order <5, which takes into account the specific geometry of such beams. A deformation v 
is split into an elementary deformation and a warping. The elementary deformation is the analog of a 
Bernoulli-Navier's displacement for linearized deformations replacing the infinitesimal rotation by a rotation 
in 50(3) in each cross section of the rod. Each part of the decomposition is estimated with respect to the 
L 2 norm of the distance from gradient v to 5*0(3). This result relies on revisiting the rigidity theorem of 
Friesecke- James-Miiller in which we estimate the constant for a bounded open set star-shaped with respect to 
a ball. Then we use the decomposition of the deformations to derive a few asymptotic geometrical behavior: 
large deformations of extensional type, inextensional deformations and linearized deformations. To illustrate 
the use of our decomposition in nonlinear elasticity, we consider a St Venant-Kirchhoff material and upon 
' various scaling on the applied forces we obtain the r-limit of the rescaled elastic energy. We first analyze 

the case of bending forces of order S 2 which leads to a nonlinear inextensional model. Smaller pure bending 
forces give the classical linearized model. A coupled extensional-bending model is obtained for a class of 
■ forces of order S 2 in traction and of order S 3 in bending. 

I. Introduction 

This paper pertains to the field of modeling the deformations of a thin structure who has a curved 
rod-like geometry with a few applications to elastic rods. Let us consider a curved rod of fixed length and 
with cross sections of small diameter of order S. Let us denote by S3 the arc length of the middle line of the 
rod, by 111(33), 112(53) two normal vectors of this line and the corresponding coordinates by s — (s\, s 2 , s 3 ). 
In this setting, the aim of this paper is twofold. In a first result, we show that a deformation v of such a rod 
can be decomposed as the sum of an elementary deformation and of a residual one as follows (see (II. 2.1)): 

(1.1) v{s) = V{s 3 ) + R(s 3 )(s 1 n 1 (s 3 ) + s 2 n 2 (s 3 )) +v(s). 

In the above decomposition, the field V(s 3 ) is the mean ofv over each section and R(s3)(sini(s3) + S2ri2(s3)) 
is the rotation of the same section, meaning that R(s3) € 5*0(3) (the special orthogonal group i.e. the set of 
orthogonal 3 x 3-matrices with determinant equal to 1). The residual field v(s) represents the warping of a 
section. The main interest of our decomposition is the fact that each term is estimated with respect to S and 
the L 2 -norm of the distance between Vv to 50(3). In order to obtain such decompositions, we first adapt the 
proof of the so called "Rigidity Theorem" established by Friesecke- James-Miiller in [11]. Our improvement 
only consists in evaluating the dependence of the quantity which measure the distance from the gradient of 
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a deformation (denned on an open set 0) to SO (3) in terms of two geometrical parameters characterizing 
(see Theorem II. 1.1). As far as thin structures are concerned, the main interest of this result is the possibility 
to slice the considered structure into small pieces for which the two geometrical parameters are uniformly 
controlled. This point is particularly helpful for a curved rod with a variable curvature which is the case 
investigated in the present paper. This allows to define the elementary deformation as a continuous field 
and to derive estimates on V, R, v and on the distance between Vv and R. These estimates first permit to 
identify a few known critical orders for the quantity ||dist(Vu, SO(3))\\l2 with respect to S (see [15], [19], 
[20]). Then, we explicitly investigate two cases namely where ||dist(Vi>, /SO(3))||x,2 is of order 8 2 and 8 K 
where k is a real number strictly greater than 2. Let us emphasize that the decomposition of v together 
with the estimates on V, R, v allow to identify the limit of the Green-Lagrange strain tensor in terms of 
the limit of the components of the decomposition of v. Moreover this decomposition of a deformation is, in 
some sense, stable with respect to the limit process with respect to S, which can be seen as a justification of 
this splitting of v. 

The second type of results concerns the asymptotic behavior of the deformations of elastic rods when 5 
goes to 0, assuming that the elastic energy is comparable to ||dist(Vw, S"0(3)) 1 and more precisely for a 
Saint Venant-Kirchhoff 's material. We consider an elastic rod submitted to dead forces (which are assumed 
to be volume forces to simplify the computations but this is not essential) . We strongly use the decomposition 
(1.1) to choose the scaling for the applied forces. In order to obtain an elastic energy of order 6 2k with k > 2, 
we are led to split the forces into two types: order 6 K ~ 1 for the loads with mean equal to over each cross 
section and order S K for general loads. We mainly investigate the cases k = 2 and k > 2. Then we also 
use our decomposition to identify the limit energy through a L-convcrgcnce argument in both cases. Let us 
briefly summarize the obtained results. 

In the case k = 2, we obtain a minimization problem which depends only on the fields V and R (and 
indeed on the forces and the boundary conditions of the 3D problem) which corresponds to the nonlinear 
energy for inextensible rods obtained in [15] and [19]. Moreover if the rod is clamped on one (and only one) 
of its extremities, we show that this minimization problem is equivalent to an integro-diffcrcntial problem 
for R and that for small enough forces there is uniqueness of the solution. 

In the case k > 2, the limit minimization problem corresponds to the standard linear bending-torsion 
energy which is also obtained in the case n = 3 in [15] and [20]. 

We also examine a situation where the forces satisfy a specific geometrical assumption (which corre- 
sponds to pure traction-compression for a straight rod) but are of order <5 K_1 (n > 3) and nevertheless which 
leads to an elastic energy of order S 2k . We obtain a linear limit model for cxtcnsional displacement in the 
elastic ID rod (with an elastic limit energy already derived in the case of a straight rod and a 3D energy of 
order <5 6 in [15] and [20]). 

As a general reference on elasticity, we refer to [7] and [3] . A general introduction to the mathematical 
modeling of elastic rod models can be found in [2], [24], see also e.g. [1], [16]. For the justification of 
rods or plates models in nonlinear elasticity we refer [1], [8], [12], [14], [15], [18], [19], [20], [21], [22], [23]. 
For a general introduction of T-convergence we refer to [9]. The rigidity theorem and its applications to 
thin structures using T-convergence arguments can be found in [11], [12], [19], [20]. For the decomposition 
of the deformations in thin structures, we refer to [13], [14] and for a few applications the junctions of 
multi-structures and homogenization to [4] , [5] , [6] . 

The paper is organized as follow. Section II is devoted to introduce the decomposition (1.1) of the 
deformations in a thin curved rod and to establish the estimates on V, R, v. In Section III, after rescaling 
the rod and the various fields with respect to S, we investigate the limit of the Green-St Venant tensor in the 
two cases ||dist(Vw, SO(3))\\ L 2 ~ S 2 and ||dist(Vw, SO(3))\\ L 2 <~ 8 K for k > 2. In Section IV, we consider an 
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clastic curved rod made of a St Vcnant- KirchhofT 's material (see IV. 1.9). After rescaling the applied forces, 
we identify the limit energy (as 5 goes to 0) through a T-convergence technique. The T-limit is a functional 
of V and R if k = 2 and of the displacement field hi and of an infinitesimal rotation field 1Z A x if k > 2. 
Then, a specific choice of applied forces leads to a linear extentional model. Section V is devoted to give an 
equivalent formulation of the limit minimization problem obtained in the nonlinear case k = 2 which leads 
to a partial uniqueness result. At least an appendix at the end of the paper details a few technical points 
concerning the interpolation between two rotations and a density result. 

II. Decomposition of a deformation in a thin curved rod 

In this section, we derive a decomposition of the type 1.1 for a deformation v of a curved rod together with the 
estimates given in Theorem II. 2. 2. In order to obtain these results, we first adapt the proof of the "Theorem 
of Geometric Rigidity" established in [11]. As mentioned in the introduction we essentially evaluate the 
dependence of the quantity which measure the distance from Vv to SO (3) in terms of two geometrical 
parameters characterizing the domain. This is the object of Subsection II. 1. Then Subsection II. 2 is devoted 
to establish the estimates on the terms of the decomposition of v with respect to ||dist(Vu, SO(3))\\l2 (see 
Theorem II. 2. 2). The techniques are similar to the ones developed for small displacements in [13] and [14]. 
At least, in Subsection II. 3 where the rod is assumed to be clamped at least on one of its extremities, we 
deduce estimates of v and Vv in terms of | |dist(Vw, SO(3))\ \^ . 

II. 1. Estimating the constant in the Theorem of Geometric Rigidity 

We equip the vector space M n ofnxn matrices with the Frobenius norm defined by 



A - ( a y)i<ij<„' 



IAIN = 
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Recall that an open set f2 of R™ is said to be star-shaped with respect to a ball B(0;R\) if for any x € 
B(0;R\) and any y <G f2 the segment [x, y] in included in f2. 

Theorem II. 1.1. Let £1 be an open set o/R™ contained in the ball B(0;R) and star-shaped with respect to 
the ball B(0; Ri), (0<Ri<R). For any v e (H 1 ^))", there exist R <= SO(n) and a e R n such that 



(77.1.1) 



||V«-R|| (i2( n))»x» <C\\dist(Vv; S0(n))\\ L 2 {n) , 
\\v - a- Rx|| (L 2(n))« < CR\\dist{\7v; SO(n)) \\ L ^n), 



R 

where the constant C depends only on n and —. 

Proof of Theorem II. 1.1. The proof of the first inequality in Theorem II. 1.1 is identical to the proof of 

Theorem 3.1 in [11] if we show that the constants which appear in the three main points of this proof only 
R 

depend upon n and — . These three main arguments arc first an approximation lemma, then a specific 
Ri 

covering of Q and finally a Poincarc-Wirtinger's type inequality. In particular, we explicitly construct a 

R 

covering of f2 which can be used in the proof of Theorem 3.1 of [11] and which only depends of — and n. 

Ri 

We begin with the following lemma which just specifies the dependence of the constants in Proposition 
A.l of [11]. 

Lemma II. 1.2. Let n > 1 be an integer and 1 < p < oo be a real number. Let Q be an open set of R" 
contained in the ball B(0; R) and star-shaped with respect to the ball B(0; Ri), (0 < R\ < R). There exists 
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a constant C = C(n,p,R/Ri) such that for any function v £ (W 1 ' p (£l)) n and for any real number A > 
there exists a function w £ (W 1 '°°(Vt)) n such that 

W ||Vu;||(ioo(n))nxn < CA 

(ii) \{ xe n ; v(x)^w(x)}\<% [ IHV^iyda; 

* p J{xen;\\\vv(x)\\\>\} 



(Hi) \\Vv - Vio||(iP(n))nxn <C 



\\\Vv(x)\\\ p dx 

i€f!; |||Vw(x)|||>A} 



Proof of Lemma II.1.2. Let us denote by B n = B(0; l) the unit ball of R™ and set S n = dB n . The 
proof of Lemma II. 1.2 is given in [11] except what concerns the dependence of the constant in the inequalities 
with respect to the geometrical parameter R/ R\ which will be extensively used in the sequel. We recall the 
following result proved in [11] ( Proposition Al; see also Evans and Gariepy [10], Section 6.6.2 and 6.6.3): 
there exists a constant Co which depends on n and p such that for any function v in (W 1 ' p {B n )) n and for 
any real number A > there exists a function w £ (W 1 ' oc (B n )) n such that 



(77.1.2) 



' ll^ 7 y w ll(L~(B„))"X" < C A 
{y€B n ; v(y)^w(y)}\<^ [ _ |||V^(y)|||^ 

A p Jfv9B„ ■ [IIV„uMII[>A* 



\\V y v - Wyw\\ {LP{Bri))nxn < C / _ \\\Vyv(y)\\\ p dy 

^ J{y^B n ;\\\Vyv(y)\\\>X} 



Since is, in particular, star-shaped with respect to the origin O, for any direction s of S n the ray issued 
from O and with direction s meets the boundary dfl on a unique point P(s). In order to transform the ball 
B n into the set fi, we first introduce the function F from S n into R + by 

Vs£5„, F(s) = \\OP(s)\\ 2 , 

where ||.||2 denotes the euclidian norm on R™ 
Now the function G from 1" into R™ is defined by 

G : y£R" — -> I MMM 

[O ]£y = 0. 

This function G is one to one from R™ onto R™ and maps B n onto 0. Moreover, due to the geometrical 
assumptions on Q, the function G is Lipschitz-continuous and satisfies the following inequalities for almost 
any y £ R" 

(77.1.3) J2iCi<|||V y G(y)|||<i2C 2 , § < ^ (x)\\\ < < | det (V w G(y)) | < R n C 2 

where the constants C\ and C 2 depend on n and R/R\. The proof of the above estimates is left to the reader 
(see also [13] and [14]). 

Let v £ {W 1,p (£l)) n . We define the function v — vo G which belongs to (W l p (B n )) n and we have for almost 
any y £ B n 

Vyv(y) = V x v(G(y))V y G(y). 
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Taking into account (II. 1.3), we obtain 



(^■1-4) fliCslllV^GO/))!!! < |||V y t;(y)||| < JK7 4 |||V x t;(G(y))||| 

where C3 and C4 depend on n and R/R\. Using the estimates on the jacobian given by (II. 1.3), we deduce 
that 

!R P R p 
<?5-^||V^||^ p(a))nxn < ||V w «||^ p(Bb))bxb < C 6 — \\V x v\\ p {LP{Q))nxn for 1 < p < 00 
G5-Rl||V x t>||(Loo(Q))nxn < ||Vj?)||(£oo( Bti ))nxn < C 6 R \ | V K V | | (£oo (O) ) n X n 

where C5 and Cg depend on n and R/R±. Now we apply the result recalled at the beginning of the proof 
so that for any A > 0, setting A = C^RX, there exists a function w £ {W 1 - aa (B n )) n such that (II. 1.2) holds 
true. Let us set w — w o G^ 1 which belongs to (W /1 ' 00 (rJ))™. Thanks to (II. 1.2) and (II. 1.5) we have 

IIV7 n ^ C0C4 R 

||V,HI(L W x„< — = — -A, 

and i) is proved. We use (II. 1.3) and (II. 1.4) to obtain 



\{x£fl; v(x)^w(x)}\ <C 2 R n \{yeB n ; v(y) ^ w(y)}\ 

/ . . IIIV^IH^ < I |||V., r ( ,,|||<V/, 



and m) is established. Now we prove Hi). We have for A and w satisfying i) and ii) 

f \\\V x v(x) - V x w(x)\\\ p dx = f \\\V x v{x) - W x w(x)\\\ p dx < 2 P f {\\\V x v{x)\\\ p + |||V x t«(a;)||| p }da; 

<2 P [ {\ p dx+ \\\V x w(x)\\\ p }dx + 2 p [ \\\V x v(x)\\\ p dx 

< C [ X p dx + 2 p [ 1 1 1 V x v(x) 1 1 \ p dx <C f 1 1 1 V x «(a;)| | \ p dx 

Jvjtw J\\\V x v(x)\\\»> J\\\V x v(x)\\\»> 

Finally we obtain 

||V x t> - V a; w|| (LP( o))«x n < C / |||V x v(x)||| p <te, 

J{xeO;|||V a ,i;(x)|||>A} 

where the constant depends on n, p and R/R\. This concludes the proof of Lemma II. 1.2. □ 

We now turn to the second argument in the proof of Theorem II. 1.1, namely the specific covering of ft. 
In the following we construct a covering Q of with cubes of the type Q(a, r) = 0+] — r, r[ n , r > satisfying 
the following properties: 

* for every Q(a, r) £ Q, the cube Q(a, 2r) is included in f2, 

R 

* there exists a constant C(n, — ) such that 

Ri 

(77.1.6) VxeQ(a,r)eQ, r < dist^x^Cl) < C(n, ^-)r, 

Ri 



* there exists a finite integer N(n) such that for every cube Q(a,r) of the covering Q, the number of 

cubes of the type Q(b,2r') which meet Q(a,r), where Q(b,r') belongs to Q, is at most N(n). 

R 

For i € N and k € N, we set n = — — = and 

TZk = | a € r fc Z" | a = r fe (zi, i 2 , • . • , i„), i p £Z and oddj. 
The covering Q of is constructed by induction as follows : 

* consider all the cubes Q(a,r ) 7 a <G TZo, such that Q(a, 2r ) C and denote by Qo the family of these 
cubes Q(a,ro) and by Wo = [J Q(a,r ), 

Q(a,ro)GQo 

* in step k > 1, consider the cubes Q(a, rfe), a e 7£fc, such that Q(a, r^) C 0\MoU...U Wfc_i and such that 
Q(a, 2rfe) C SI, and denote by Qk the family of these cubes Q(a, rfe) and by Uk = [J Q(a, rfe). 

Q(a,r k )eQu 

We denote by Q the countable family of all the cubes constructed through the above process. 

The above explicit construction permits to show that the covering Q verifies the required properties (as an 

R R 

example we can take C(n, — ) = b\/n— and N(n) = 2™ +3 ). 

ill Rl 

As far as the third argument in the proof of Theorem 1.2.1 is concerned, we now recall the following 

Poincare-Wirtinger's inequality (see [13] for a proof and various applications). Since SI is contained in the 

ball i?(0; R) and is star-shaped with respect to the ball B(0] Ri), there exists a constant C which depends 
R 

on n and — such that for any </> G _ff 1 (Sl) (see [13]) 
Ri 

(H.1.7) U - M{4>)\\ LHn) < C\\ P V4>\\ {L 2 {n)) n 

where M.{4>) is the mean of </> over SI and p(x) — dist(x, dQ). Using Lemma II. 1.2, the specific covering of il 

described above and the Poincare-Wirtinger's inequality (II. 1.7) permit to reproduce the proof of Theorem 

3.1 in [11] in order to obain the first estimate of our Theorem II. 1.1 with a constant which depends only on 
R 

n and — . To end the proof, we apply inequality (II. 1.7) to the field v(x) — Rx and we use the first estimate 
Ri 

in Theorem II. 1.1. □ 

II. 2. Decomposition of the deformation in a curved rod. Estimates 
II. 2.1. The geometry 

Let us introduce a few notations and definitions concerning the geometry of a curved rod (see [13], [14] 
for a detailed presentation). 

Let ( be a curve in the euclidian space R 3 parametrized by its arc length S3. The current point of the 
curve is denoted M(s 3 ). 

We suppose that the mapping S3 — > M(ss) belongs to (C 2 ([0, L})) 3 and that it is one to one. We have 

dM -t lltll -1 

1 — - F 2 - 1, 

ds 3 

where || • H2 is the euclidian norm in M 3 . 

Let ni be a function belonging to (C 1 ([0,L])) 3 and such that 

Va 3 e[0,L], ||ni(s 3 )|| 2 = 1 and t(s 3 ) • rn (s 3 ) = 0. 
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We set 

n 2 = t A ni. 

In the sequel, oj denotes a bounded domain in M 2 with lipschitzian boundary (while obviously, ZJ denotes 
the closure of oj). We choose the origin O of coordinates at the center of mass of w and we choose the coor- 
dinates axes (0;ei) and (0;e 2 ) as the principal axes of inertia of w, so that / xidxidx 2 — \ x 2 dx\dx 2 = 

J U J LJ 

I x\x 2 dx\dx 2 — 0. The reference cross-section lo$ of the rod is obtained by transforming u with a dilatation 

J u 

of ratio S > and we set 

Vlg — u$ x (0, L). 

Introduce now the mapping $ : M 2 x [0, L] — > M 3 defined by 

$ : {si,s 2 ,s 3 )\ — > M(s 3 ) + s 1 n 1 (s 3 ) + s 2 n 2 (s 3 ) 

There exists So > depending only on £, such that the restriction of $ to the compact set Q$ is a C 1 — 
diffcomorphism between Qg and $(il<5). Moreover, there exists two positive constants c and ci such that 

V<ye(0,5 ], VseUs, c< |||V$(s)||| < ci. 

Definition II.2.1. For S e (0,<5 ], the curved rod Vs is defined by 

v s = $(n s ). 

The cross-section of the curved rod is isometric to oos- In Vs, the point M(s 3 ) is the center of gravity of 
the cross-section &(ujs x {S3}) and the axes of direction ni(s3) and n 2 (s3) are the principal axes of this 
cross-section. 

Notation. Reference domains and running points. We denote x and s respectively the running point of Vs 
and of so that x = <&(s). 

A deformation v defined on Vs can be also considered as a deformation defined on fig which we will also 
denote by v, as a convention. As far as the gradients of v are concerned we have V s v = V^v-V^ 

VS e (0,5 ], c|||V^(x)||| < |||V S «( S )||| < q||V x w(a:)||| 

where the constants are positive and do not depend on S. 
II.2.2. The elementary deformation 

In this subsection, we show that any deformation v € (H 1 (Vs)) 3 of the rod Vs can be decomposed as 
(using the above convention) 

(77.2.1) v(s) = V(s 3 ) +R(s 3 )(s 1 n 1 (s 3 ) + s 2 n 2 (s 3 )) +v(s), setts, 

where V belongs to (i7 1 (0,i)) 3 , R belongs to (77^0, 7)) 3x3 and satisfies for any s 3 e [0,7]: R(s 3 ) e 5*0(3) 
and v belongs to (H^Vs)) (or (H^flg)) using again the same convention as for v). Let us give a few 
comments on the above decomposition. The term V gives the deformation of the center line of the rod and it 
is indeed a function of the arc length S3. The second term R(s 3 ) (sini(s 3 ) + s 2 n 2 (s3)) describes the rotation 
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of the cross section (of the curved rod) which contains the point M(ss). The sum of the two first terms 
V(s3) + R(s3)(sini(s3) + s 2 n 2 (s3)) is called an elementary deformation of the rod. 

II. 2. 3. The main theorem 

The following theorem gives a decomposition (II. 2.1) of a deformation and estimates on the terms of 
this decomposition. 

Theorem II. 2. 2. Let v e (i? 1 ^)) 3 , there exists an elementary deformation V + R(sini + s 2 n 2 ) and a 
warping v satisfying (II. 2.1) and such that 



(77.2.2) 



Ml(z,2(n 5 ))3 < C5\\dist(T7 x v,SO(3))\\ L 2(p s) 
|V s zJ||( L 2( J2iS ))3x3 < C\\dist{V x v, SO(3))\\ L 2( r5) 



dn 



ds 3 



< 



C 



(L 2 (0,L)) 3 > 



-\\dist(V x v,SO(3))\\ L 2 (Vs) 



ds 3 

\V x v-R\ 



C, 



(ij()in3 < T l|^(V^50(3))|| L2(n) 



(!!s)) 3x 3 < C\\dist{V x v,SO{i))\\ L 2 {Vs) 



where the constant C does not depend on S. 

\2L Li L 
Proof. Let N be an integer belonging to — -, — and let < a < L . 

L 3 

We have 5 < — < -5. Let R > 1 be such that the reference cross-section uj is contained in the ball B(0; R). 

Then the domain Cls.a = wjx]a,a+ — [ has a diameter less than 3RS. In the sequel we will work with the 
portions Vs, a of the rod Vs defined by 

P S , a = *(fi 4 ,a). 

As in [13] we distinguish two cases. 

First case. If uo is star-shaped with respect to a ball of radius R\ < 1/2, it is shown in [13] that each Vs. a 

R S 

is star-shaped with respect to a ball of radius — — and we are in a position to apply Theorem II. 1.1 to the 

8 

R 

function v into each part Vs, a for which the ratio — is independent of 5. As a consequence, there exist 

R\ 

R Q e SO(3) and a Q e K 3 such that 



(77.2.3) 



\V x v - R Q ||( L 2 (n>c<)) 3x3 < C||dist(V x u;S'0(3))|| z ,2( 7 , i5a) 

\v-&a- R„ (x - M(a))|| (L 2 (na)) 3 < C5\\dist(V x v;SO(3))\\ L 2(p 6ia y 



The constant C does not depend on a and S. 

Second Case. Let us consider the general case, where the cross-section is a bounded domain in R 2 with 
lipschitzian boundary. There exists a finite sequence of open sets . . . ,u}^ such that 



CO 



- u « ( °. 



ojs= (J u\\ 



KKK 



1<1<K 



and such that every oj^ is star-shaped with respect to a disc of radius Ri, < R\ < 1/2. Moreover, the 
open set uj is connected, then there exists 7? 2 s]0, Ri] such that uj^ n uj^ contains a disc of radius 7? 2 if 
the intersection is not empty. 
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The domain fl^ a — cj^x]a, a + — [ is star-shaped with respect to a ball of radius RiS. As in the first 



case, there exist R„' e SO (3) and al, 1 ' e M 3 such that 



(0 



|Vx« - Ri ll (i2(P (o )) 3X3 < C||dist(V xU ; 50(3))|| L2(T , (! , . 



| v -ag)-R(|) (s -M(a)) || 



(i 2 (Pi'i)) 3 



< C5||rfist(V x w;6'0(3))|| 



The constant C does not depend on a , 5 and /. 

If a/ r ' n ix>( s ) 7^ the portion pjfl O contains a ball of radius R26/8. This allows us to compare the 
elementary deformations a„' + [x — M{aj) and a«' + R«' (x — M(a)) in this ball. We obtain 



C 



RL r) -Ri s) ||| 2 < ^\\dist(y x v;SO(3))\ 



\ a ^~ a ^\\ 2 < j\\d lS t(V x v;SO(3))\ 



where the constant only depends on R , i?i and i? 2 - 

Setting R Q = R„' and & a = ai"^ and proceeding step by step with respect to I, we finally deduce that 
(II. 2. 3) holds true with a constant C which does not depend on a and 5 as in the first case. 

Now we consider two splittings of Vs by considering two sets of arc length 

L L 
ak = k N> k = >---> N > Pk^a k + — 1 k = 0,...,N-l. 

We consider the elementary deformations a Qfc + R Qfc (x — M(otk)) and a.p k + R^ fc (x — M((3k)) of the 
portions Vs,a k an d V&fik which satisfies estimates (II. 2. 3) with a constant independent of k. Considering 
Vs,a k n Vs,p k and Vs,a k+1 H T 5 ^, we can compare R ak and R« fc+1 . We obtain 



(77.2.4) 



iV-l 

V- 

^ iV 



7/7V 



2 <^\\di S t(V x v;SO(3))\\l 2(Vs) , 



where R QJV = R aw l and where the constant C is independent of 5. 

Now we are in a position to define the elementary deformation associated to v. We set for S3 € [0,7] 



(77.2.5) 



V(s 3 ) = — ^-r / u(si,S2,S3)dsids 2 - 



In order to define R we use the following argument whose proof is postponed to the appendix. There 
exists a field of matrices R belonging to (i7 1 (0, 7)) 3x3 , with R(s 3 ) € 50(3) for all s 3 e [0,7], such that 
R(afc) = R Qfc for k = 0, . . . , N and 



(77.2.6) 

Indeed the field v is defined by 



dR 



ds 3 



< 



(L 2 (0,L)) 3 >< 3 



N-l 

k=0 



N 



R"fc+i R"fc 



L/N 



(77.2.7) 



v(s) = v(s) - V(s 3 ) - R(s 3 )(sini + s 2 n 2 ) for a. e. s e O5. 



The third estimate of (II.2.2) follows directly from (II.2.4) and (II.2.6). 



9 



From (II.2.3) we obtain 

JV-l 

l|VxW-R a J|^ (7 , 5aJ) 3x3 <C||dist(V x «;SO(3))||i a(7 , 4) 

fc=0 
N-l 

\\v-a ak -R ak (x~M(a k ))\\f L2(VeaJ)3 < C5 2 \\dist(\7 x v; SO(3))||| 2(n) . 

fc=0 

Now, we take the mean value over the cross-sections of Cts, then using the definition of V we deduce 



(77.2.8) 



N-l 



(77.2.9) ^||V-a Qfc -R Qfc (M-MK))|| (i2(QfcQfc+i))3 <C||^(V^;50(3))||i 2( ^ 



fc=0 



Thanks to the definition of R and the third estimate in (II. 2. 2) we get 
(77.2.10) 



N-l 



c 

ll R _ R «* H(L*(a»,a» +1 ))3*s < 72 \\dist(V x V, 50(3)) | || 2(n) . 
fe=0 



Consequently (II. 2. 9) and (II. 2. 10) give the first estimate in (II. 2. 2) while (II. 2. 10) leads to the last estimate 
in (II.2.2). 

Due to the definition of $, we have — — = n Q and — — = t + si^— !- + S2-r— ■ so that the relation 

ds a asz ds 3 ds 3 

V s v = V x v.V& leads to 



. -i -i \ dv „ „ „ / dri! dn 2 \ 

77.2.11 7-=V I »n 1 — = V K wn 2 — = t + Sl — + s 2 — . 

osi 0S2 0S3 V as3 as3 / 

Then inserting (II. 2. 10) into (II. 2. 8) gives in particular 
(77.2.12) 



dv 



dv 



R(t + si- h s 2 -^— 

OS3 ds3 ds 3 



(i 2 (^)) 3 



< C\\dist(V x v;SO(3))\\ 2 L2{Vi 



Integrating first over w<5 x {53} in (II. 2. 12) leads to the forth estimate of (II. 2. 2) (recall that / S\dsids2 

J u 

I s 2 dsids2 = 0). It remains to show the estimate on V S U. From (II. 2. 11) we get 

J u 



(77.2.13) 

and then with (II.2.8) and (II.2.10) 
(77.2.14) 

dv 



dv dv ._ 

— - = - R) ni -— = {V x v - R) n 2 
asi as 2 



dv 


2 


dv 


2 


dsi 


+ 


dS2 


(i 2 (^)) 3 



< C\\dist(V x v;SO(3))\\ 2 L2{Ps) . 



Now we estimate 



<9s 3 



(L2(n a ))3 



We have from (II.2.7) and (II.2.11) 



77.2.15 — - (V x w - R) (t + si— i + s 2 — i) - — «mi + s 2 n 2 ) - — + R t. 

as 3 v 7 v ds 3 ds 3 ' ds 3 v ' as 3 
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Proceeding as above to bound the first term and using the third, the forth and the last estimates of (II. 2. 2) 

dv 

to control the last two terms of the above inequality permit to obtain the estimate on — — given in (II. 2. 2). 

ds 3 

□ 

In order to split the bending and the stretching of the rod, we now introduce the following quantities: 

(77.2.16) Vs 3 G [0, L], V B (s 3 ) = V(0) + / * H(z)t(z)dz, V s (s 3 ) = V(* 3 ) - V B (s 3 ). 

Jo 

Let us notice that Vb is the bending deformation of the middle line while Vs is the stretching deformation. 
Due to the third estimate in (II. 2. 2) we have 

(^■2-17) \\Vs\\ { m { o,L))s ^ j\\dist(V x v,SO(3))\\Li(Vs)- 

Inserting the definition (II. 2. 16) into the decomposition (II. 2.1) gives 

(77.2.18) v(a) = V B (s 3 ) + R(s 3 )(s ini (s 3 ) + s 2 n 2 (s 3 )) + V s (s 3 ) + v(s), selij. 

Estimates (II. 2. 2) and (II. 2. 17) permit to interpret the part Vs(s 3 ) + R(s 3 )(sini(s 3 ) + s 2 n 2 (s 3 )) of 
the decomposition (II. 2. 18) as an approximation of the parametrization of the deformed rod at least if 
||dist(V^,50(3))|| i2(n) «8. 

II.2.4. The boundary condition 

Let us denote by I 3 the unit 3x3 matrix and by Id the identity map of R 3 . 

In this subsection, we derive boundary conditions on the terms of the elementary deformation given 
by Theorem II. 2. 2. Indeed these conditions depend on the boundary conditions on the field v. We discuss 
essentially the usual case of a clamped condition on the extremity of the rod defined by 



r , 5 = $(UJ6 x {0}). 

Then we assume that 

v(x) = x on ro,5. 

In the following we show that the elementary deformation V(s 3 ) + R(s 3 )(sini + s 2 n 2 ) given by Theorem 
II. 2. 2 can be chosen such that v = on the corresponding boundary ivs x {0}. Due to the definition (II. 2. 5) 
of V, we first have 

V(0) = M(0), 

(recall that the point M(0) is the beginning of the middle line of the rod Vs', see the notations in Subsection 
3.1). Now we prove that in Theorem II. 2. 2 the choice R(0) = I 3 is licit as a boundary condition for the 
matrix R(s 3 ). Estimates (II. 2. 3) for the first portion Vs, ao imply 



K, ., 0) - a Q1 - R Q0 (s ini (0) + s 2 n 2 (0)) | \ 2 {L 2 M)3 < C5dist(V x v: SO(3)) ||| 2(n o 



Using now the boundary condition written in the equivalent form v(si, s 2 , 0) = M(0) + Sini(0) + s 2 n 2 (0) in 
the above estimate and using again / s\ds\ds2 = / s 2 dsirfs 2 = / sis 2 dsirfs 2 = lead to 
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IIIRao -I3III 2 < ^\\dist(V x v;SO(3))\\ 2 L2{Psao y 

As a consequence we can substitute R Qo by I 3 in the construction of the function R(s 3 ) without altering 
estimates (II. 2. 2) so that R(0) = I3. Indeed this leads to v = on the boundary L05 x {0}. The above 
arguments can be easily adapted if the imposed deformation on Toj is of the form v{x) = A + P(x — M(0)) 
where A is the deformation of the point M(0) and P is 3 x 3 matrix. This leads to two boundary conditions 
of the type V(0) = A and R(0) = Q where, as an example, the rotation Q minimizes the distance from 
P to 50(3). In the same spirit, if the rod is clamped on the two extremities of Vs, one can modify the 
construction of the function R(sa) in such a way that R(0) = R(£) = I3 keeping (II. 2. 2) true. 

II.3. 77 1 - Estimates 

Throughout the paper we now assume that the boundary Tq^s is clamped so that 

(77.3.1) V(0) =M(0) = Vs(0), V s (0) = and R(0) = I3 

and we denote by C a generic constant independent of S. 

We derive a first 77 1 - estimates using directly (II. 2. 2) and the fact that ||Rt||2 = 1- It gives 

(77.3.2) ||£|| (i2(oi))3 <c(l + i||dist(V^50(3))|| L2(n) ). 
Using the boundary condition (II.3.1), it leads to 

(^■3.3) l|V|| (L2(0 , L)) 3 <c[l + - s \\di S t(V x v,SO('S))\\ LHVs) ). 

Inserting (II. 3. 2), (II. 3. 3) into (II. 2.1) and using the estimates of (II. 2. 2), we deduce that 



(77.3.4) 



\M(L2(V 5 )) 3 + \\VxV\\(L2(Vs)) 3 * 3 < C ( S + ||dist(Va:U, SO(3))\\ L 2(p 5 ^ , 

J|«-V|| (i2(n)) 3<CT(<5 + ||dist(V^,50(3))|| i2(n) ). 



Notice that using (II. 2. 2) also leads to the following estimates 

C 



R I 



(77.3.5) 




3 ll(L2(o,D)3x3 ^ ^lldist^w, SO(3))\\ L 2 (Vs) 



(L2(oi))3 <g||dist(V^,50(3))|| i2(n) . 



Since t = 4— the last inequality of (II. 3. 5) together with the boundary condition (II. 3.1) give 

(^•3-6) \\V - M\\ (L2{0 L))3 < ^||dist(V x «,SO(3))|| L2(7 , 4) . 

From (II.2.1), (II.3.5), (II.3.6) and estimates (II.2.2), we deduce that 

C 

(77.3.7) \\v - 7 d || (L 2 (n))3 + \\V x v - I 3 ||( L 2 (7 , 4 ))3X3 < — ||dist(V ;E w,50(3))|| L 2 (n) . 
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From (II.2.2) and (II. 3. 5) we also have 

(^•3-8) \\R-l4 (H1{0 L))3x3 < ^||dist(V^,SO(3))|| i2(n) . 

From (II.2.2) and (II.3.8) and the fact that (v - I d ) - (V - M) = (R - I 3 )(sini + s 2 n 2 ) + v, we obtain 

(77.3.9) \\(v - I d ) - (V - M)|| (W)) 3 < C||dist(V x t>, SO(3))\\l H v s ) 

and 

||V X « + (V x vf - 2I 3 || (L2(n5)) 3x3 < q|dist(V x «, SO(3))\\ L 2 (V5) + CS\\K + R T - 2I 3 || (L2(0 i))3x3 . 

Due to (II.3.8), the continuous embedding of iJ^O.L) in C°([0,L]) and the equality R + R T 2I 3 
R T (R - I 3 ) 2 we finally obtain 



(77.3.10) { c 



V x v + (V x vf -2I 3 \\ (L2(Px))3x3 <C\\dist(V x v,SO(3))\\ L 2 (V5) 

^\\dist(V x v,SO(3))\\l H v s y 

III. Asymptotic behavior of a sequence of deformations 

In view of the first estimate in (II. 3. 5) we can distinguish three main cases for the behavior of the 
quantity ||dist(Va;U, SO(3))\\ L 2^ Vs - ) (which will be a bound from below of the elastic energy) 



||dist(Vxt>,SO(3))|| L2(7 > 4) 



0(5 K ), 1<k<2, 



k > 2. 




This hierarchy of behavior for ||dist(V x w, £0(3)) || £2^) has already be observed in terms of clastic energy 
in [15]. 

In this section we investigate the behavior of a sequence of deformations of Vs for ||dist(V x u, SO(3))\\l2(p s } = 
0(5 K ), for k > 2. Actually, the first estimate in (II. 3. 5) is useless if k < 2 then we can not analyze this case 
using the decomposition (II. 2.1). As usual we first rescale Sis in order to work over a fix domain in Subsection 
III. 1 . In Subsection III. 2 we investigate the case k = 2 while in Subsection III. 3 we deal with k > 2. Let us 
emphasize that we explicit the limit of the Green-St Venant's tensor in both cases. Subsection III. 4 gives a 
few comparisons with the linearized deformations. 

III.l. Rescaling of Sis 

We set SI = uj x (0, L) and, we rescale Sis using the operator 

(Ils<j>)(S 1 ,S 2 ,s 3 ) = <j>(dS u 5S 2 ,s 3 ) for any {S u S 2 , s 3 ) <= SI 

defined for any function <j> defined over Sis- Indeed, if <f> € L 2 (Sls) then (Hs(f>) € L 2 (Sl). The estimates (II. 2. 2) 
of v transposed over SI are 

' ||n^|| (i 2 (n))3 < C\\dist(V x v,SO(3))\\ L 2 {P5) 
|^|| ( l W < C\\di s t(V xVl SO(3))\\ L 2 {Ps) 

(777.1.1) { dllsv 

l |^||(L 2 (n))3<q|dist(V :c «,50(3))|| i2( ^ ) 

|^^||(^(n))3 < j\\dist(V x v,SO(3))\\ L 2 {V6) 
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while estimates (II. 3. 4) and (II. 3. 7) on v become 

' ||lMI(L*(fi))3 < C(l + i||dist(V^,50(3))|| L2(n) ) 
'^'"Vw <c(s + \\di S t(V x v,SO(3))\\ L 2 {Ps) ) 



(777.1.2) 



1 dSi 1 
, dUgv . 
1 95 2 1 

1 a S3 1 



(LW <C(* + ||dist(V x i;,SO(3))|| ia(7 , 4) 



l(i 2 (n))3 



< 



c(l + i||dist(V^,50(3))|| L2( ^)) 



and 



(777.1.3) 



C 



n 5 (u-7 d )|| (i 2 (a)) s < -||dist(V ;E «,50(3))|| i 2 (n) 



9n 5 («- 


Id) 






0n«(« - 


Id) 


dS 2 




dii s (v - 


Id) 


ds 3 



|(^(Q))3 < y||dist(V x w,SO(3))|| i 2 (7 , 5) 
C 

l(i 2 (n)) 3 < yl|dist(V a; w,S'0(3))|| L 2 (T , (S) 
C „ 



All the estimates given in Section II. 2. 3 over Qg can be easily transposed over f2. 

III. 2. Limit behavior for a sequence such that ||dist(V x U5, 50(3))|| L 2 (n) ~ (5 2 

Let us consider a sequence of deformations vg of (Tf^T-^)) 3 such that = 7^ on T 0y g and 

(777.2.1) ||dist(V x ^,50(3))|| i 2 (n) < CS 2 . 

Indeed different boundary conditions on vg can be considered on both the extremities of Vg (see Subsection 
II. 2. 4). We denote by Vg, R5 and vg the three terms of the decomposition of vg given by Theorem II. 2. 2 
and by V B ,s and V s ,s the two terms given by (II.2.16). The estimates (II.2.2), (II.2.17), (II.3.5), (II.3.6), 
(III. 1.1), (III. 1.3) lead to the following lemma: 

Lemma III. 2.1. There exists a subsequence still indexed by S such that 

Kg^R weakly in (77^0, 7)) 3x3 
Vg — >V strongly in (^(O.L)) 3 
V B .s — > V strongly in (77^0, L)) 3 



(777.2.2) 



■jV s ,s -» V s weakly in (H\0, L))'' 



1 



■jpHgVg — 1 v weakly in (7 2 (0, 7; i7 1 (w)))' 



Moreover R(ss) belongs to 5*0(3) for any S3 e [0,7], V € (77 2 (0,7)) and they satisfy 

dV 



(777.2.3) V(0) = M(0), R(0) = I 3 , V s (0) = 0, 

Furthermore, we also have 



anc 



ds 3 



= Rt. 



(777.2.4) 



V strongly in (T/ 1 ^)) 3 , 



na(V x Va) — >■ R strongly in (7 2 (f2)) 3x3 . 
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The fields V and R which arc defined in Lemma III. 2.1 describe the deformation of the limit ID curved 
rod as a deformation of the middle line V and a rotation of each cross section R. The convergences (III. 2. 4) 
show that (V, R) is the limit of the deformation of the (rescaled) 3D curved rod and represents the elementary 
deformation of this rod. 

dV 

The last relation in (III. 2. 3) shows that - — = 1 and then the variable S3 remains the arc length of 
the middle line of the deformed configuration. As a consequence, in the present case, there is no extension- 
compression of order 1 the rod. Then the limit behavior is essentially a bending model. At least, the forth 
convergence (III. 2. 2) means that the quantity Vs which describes the stretching is of order S. 

The following corollary gives a corrector result. 

Corollary III. 2. 2. For the same subsequence as in Lemma III. 2.1, we have 



(//J.2.5) 



' 1 d — 

-^(U s (V x Vs) - R«5)n Q weakly in (7 2 (ft)) 3 , 

7(n 4 (V x « 4 )-R*)t^^(Smi + S 2 n 2 ) + ^- weakly tn (7 2 (f!)) 3 . 
ds3 ds3 



Proof. The first convergence in (III. 2. 5) is a direct consequence of (II. 2. 13) and (III. 2. 2). In order to 
obtain the second convergence, remark first that, thanks to estimates (III. 1.1) and (III. 2.1) the sequence 
-^Hsvs is bounded in i? 1 (fi). Due to (III. 2. 2), its weak limit must be equal to 0. Using now (II. 2. 15) and 
the convergences (III. 2. 2) leads to the result. □ 

To end this section, let us notice that the strong convergences in (III. 2. 4) together with the relation 
(ViWj) T ViU5 — I3 = iV x vs — Rs) T V x vg + (R5) T (V X V5 — Rs) permit to obtain the limit of the Green- 
St Venant's tensor in the rescaled domain Q 

(III. 2.6) ^U s ((V x v s ) T V x v s -I 3 ) weakly in (L 1 (ft)) 3x3 , 



where 



{II 1. 2.7) 



„ If, . . / dv . dv . dR,„ „ , dVsV 1 ' 

= 2^ 1 n2 1 1 dS 2 1 ds-^ + ^ + *f ) 



_ T ( dv . dv . dR , „ . dVs \ , , T i 

+R (m 1 w 2 1 dsi ^ + + *f ) (ni 1 n2 1 1} I 



\dS± dS 2 ds3 ds 3 



and where (a | b | c) denotes the 3x3 matrix with columns a, b and c. 

dR 

Setting w = R T v and using the fact that the matrix R T - — is antisymmetric, we can write E as 

ds 3 

(7/7.2.8) E = (m I n 2 1 1) E (m | n 2 | t) T , 

where the symmetric matrix E is defined by 



(7/7.2.9) E 



I dw lr dw dw 1 1 f dw „ dR ^ dVs ~ I 

(ds;- ni 2fe' n2 + ^' ni ) 2\ds;- t - S2 dr 3 ni - Kn2 + ^- Ilni \ 

dw 1 f dw „ dR _ dVs ^ I 

ds- 2 - n2 2{w 2 - t + Sl ds- 3 ni - nn2 + ^;- Rn2 \ 

„ dR „ „ dR „ dV s „ 
* * St— t • Rni - S 2 — t • Rn 2 + — — • Rt 

ds 3 ds 3 ds 3 
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Corollary III. 2. 3. Assume that 



V<5 e]0,6 ], dct (Vv g (x)) > for a.e. xe Vs 
then, for the same subsequence as in Lemma III. 2.1, we have 

(7/7.2.10) ||E|| (i2( n))3x3 = ||E|| (L 2 (n))3x3 < lim -^||dw*(V*t;«, SO{Z))\\ L 2 {Vs) . 

(5^0 ° 



Proof. The map A — > \J A T A is continuous from the space of the 3x3 matrices into the set of all 
symmetric matrices and we have |||\/A' r A||| = |||A|||, where ||| • ||| is the Frobenius norm. Then, the second 
strong convergence in (III. 2. 4) gives 



n^yVzt^V^) — -> I 3 strongly in (L 2 (Q)) 



3x3 



Estimate (III. 2.1) implies that the sequence —H$ (V x vs) T \7 x vs — I3J is bounded in (7 2 (£l)) 3x3 . The 
identity (V ' x v s ) T V ' x v s - I 3 = (^J {V x v s ) T V x v$ - I 3 ) ( (V x v&) T V x vg + I 3 ), the weak convergence (III. 2. 6) 
and the above strong convergence give 

^U s y(V x v s ) T V x v s - I 3 ) ^ E weakly in (7 2 (0)) 3x3 . 

We recall that for any 3x3 matrix A such that dct(A) > 0, we have dist(A, 50(3)) = |||V.4 T A - I 3 |||. By 
weak lower semi-continuity of the norm, we obtain the result. □ 

III. 3. Limit behavior for a sequence such that ||dist(V a; u 5 , SO(3))\\ L 2 {Vs) ~ 6 K for k > 2. 
Let us consider a sequence of deformations vs of (i7 1 (7 , 5)) 3 such that vs = Id on T 0y s and 

\\dist(W x vs,SO(3))\\ L 2 {P5) <CS K . 

The estimates (II. 3. 7) and (III. 1.3) lead to the following convergences: 



(777.3.1) 



Ra^Is strongly in (H\0, 7)) 3x3 , 
U s vs —> M strongly in (/^(ft)) 3 , 
[ILj(V x O^I 3 strongly in (7 2 (0)) 3x3 . 



We now study the asymptotic behavior of the sequence of displacements 

u g =v s - I d . 

Due to decomposition (II. 2.1) we write 

(777.3.2) u 6 (s) =Us{s 3 ) + (R 5 - I 3 )(s 3 )( Sini (s 3 ) + s 2 n 2 (s 3 )) +v s (s), s e Sl s , 

where %(s 3 ) = Va(s 3 ) — M(s 3 ) = (Vs,i(s 3 ) — M(s 3 )) + Vs,<5(s 3 ) and we have the following Lemma: 
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Lemma III. 3.1. There exists a subsequence still indexed by S such that 



(7/7.3.3) 



F^(R<5 - I 3 ) A weakly in (Tf^O,^)) 



3x3 



-^Z^Us — >U strongly in (i7 1 (0,7))" 



— Vs,8^V s weakly in (7^(0, 7))" 
—Usvs — 1 v weakly in (7 2 (0, L; H (co)))' 



The function U belongs to (i7 2 (0,i)) 3 , for any S3 € [0,7] f/ie matrix A(s^) is antisymmetr 
following relations hold true: 



and the 



(777.3.4) 

Moreover we have 



U(0) = V s (0) = 0, A(0) = and ^ = At, 

as 3 



(777.3.5) 



-^glLjtij — > W strongly in (77 (fi)) , 
— ^Il^Vatia) — >■ A strongly in (i 2 (£l)) 3x3 . 



Proof. The convergences (III. 3. 3), (III.3.5) and the relations (III. 3.4) follow directly from estimates (II. 2. 2), 
(II. 3. 6), (II. 3. 7) and (III. 1.3). It remains to prove that A(s3) is antisymmetric. Using the first convergence 
in (III. 3.1) and the first convergence in (III. 3. 3) we get 



— Rj(R, - I 3 ) — ► A strongly in (L 2 (0, L)) 



3x3 



The matrix R5 belongs to 50(3), hence Rj(Ra — I 3 ) = — (R5 — I3) T . Then, from the first convergence in 
(III. 3. 3), we deduce that the matrix A(s 3 ) is antisymmetric. □ 

Since A is antisymmetric, there exists a field 1Z e (i7 1 (0, L)) 3 (with TZ(0) = 0) such that for all x e M 3 



(777.3.6) 



Ax — TZ Ax. 



From (III. 3. 4) and the above equality we obtain 



(777.3.7) 



dU 

ds 3 



n At. 



The relation (III. 3. 6) means that, at the order 5 K 2 , the cross sections of V& are submitted to small rotations 
and (III. 3. 7) shows that the limit displacement is of Bernoulli-Navier's type. 

Corollary III. 3. 2. For the same subsequence as in Lemma III. 3.1, we have 



(777.3i 



^^(UsCVxVs) - Rs)n a weakly in (7 2 (fi)) 3 , 



dS a 

^ T (n 4 (V I t> 4 )-R*)t^^A(5ini+5 2 n2) + ^ weakly in (7 2 (ft)) 3 . 
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Proof. The proof of Corollary III.3.2 is similar to that of Corollary III. 2. 2, but using now (II. 2. 13), 
(II. 2. 15) and the convergences of Lemma III. 3.1. □ 

From Lemma III. 3.1 and Corollary III. 3. 2 we deduce the limit of the Green-St Venant's tensor in the rescaled 
domain fi 

(7/7.3.9) _L_n 5 ((V x ^) T V xU5 -I 3 ) weakly in (L 1 (Q)) 3 * 3 , 

where the symmetric matrix E is defined by 

„ 1 r / dv . dv . dTZ , „ dVs \ T , 

/ i i .it I &v , dv , dTZ , , dVs \ l 

+(ni|n2|t) fe'^'^^^+^ + ^r)} 

We can write E in the form 

E = (m | n 2 | t)E(m | n 2 | t) T 
where the symmetric matrix E is defined by 



(7/7.3.10) E 



/ dv 1 r dv dv ) 1 r dv dTZ dVs 1 \ 

dv 1 r dv „ dTZ dV s 1 

^• n2 2m- t+Si d7- t+ i^- n2 \ 



„ dTZ „ dTZ dV s 

\ * * -Si-— ■ n 2 + S 2 -r~ ■ ni + — — t / 

\ as 3 ds 3 ds 3 / 



From Lemma III. 3.1 and the above convergences, we deduce the analog of Corollary II. 2. 3. 
Corollary III.3.3. Assume that 

V<5 e]0,<5 ], det (Vv s (x)) > for a.e. xeV s 

then, for the same subsequence as in Lemma 111.3.1, we have 

(777.3.11) ||E|| (i 2 (0)) 3x3 = ||E|| (i 2 (0))3 x3 < }m-^\\dist(S7 x v s ,SO(3))\\ L 2 (Vs) . 



III.4. Comparison with linearized deformations 

In this subsection, we always consider a sequence of deformations v$ of (//^'Pa)) 3 satisfying v$ = Id on 
T .s and 

\\dist(V x vs,SO(3))\\ L 2 {Ve) <C5 K . 

We recall the decomposition (III. 3. 2) of the displacement us = v$ — Id- 

Let us notice that (II. 3. 5) shows that for n > 2, both the displacement and its gradient are small (with 
respect to S). One can then address the problem of comparing the limit displacement Us and the limit 
displacement in the framework of small deformations. To this end let us first recall the decomposition of 
displacement for small deformations. 

We define the strain semi- norm | • \s by setting 



1 



Vw e {H\V 5 )Y \w\ £ = \\-(V x w + (V x w) 



2\ 1 ' \ ) ||(L2(-p 4 ))3x3 



18 



Now, using the results obtained in [13], we decompose ug in the sum of an elementary displacement and a 
warping 

ug(s) = U et s(s) +us(s) = U s (s 3 ) + 7lg(s 3 ) A (sini(s 3 ) + s 2 n 2 (s 3 )) + ug(s) for a.e. s E Clg. 
The warping us satisfies the following equalities 

/ us(si,s 2 ,s 3 )dsids 2 = / us(si,s 2 ,s 3 ) A (sini(s 3 ) + s 2 n 2 (s 3 ))dsids 2 = for a.e. s 3 e (0, L). 

JU>5 J LOS 

Notice that the first term Us of the elementary displacement U e ,g is the mean value of us over the cross-section 
<&(uos x {S3}) and then is the same as in (III. 3. 2). Theorem 2.1 in [13] gives 



(777.4.1) 



|V s U5||(i,2(o i5 ))3x3 < C\us\s \\us\\(L 2 (n s )) 3 < CS\us\s 



dTZs 



ds 3 



US\£ 



<C ... 

(L 2 (0,L)) 3 S 2 



dU 5 
ds 3 



lis At 



US\£ 



<C . 

(L 2 (0,L)) 3 S 



We recall (see [14]) the definitions of the inextensional displacements and extensional displacements sets of 
the middle-line of the curved rod. 



(777.4.2) 



D In = {lle (H\0, L)f I ^ • t = 0, U(0) = o} 



D Ex = {U€ (77^0, 7)) 3 I — -ni 



ds 3 

dU dU 1 



An element of 7)/„ is an inextensional displacement while an element of D Ex is an extensional one. We recall 
(see [14]) that Us can be written as the sum of an inextensional displacement and an extensional one 

(777.4.3) Us = U I .s + U E ,8 Ui, s € D In , U E ,s e D Ex , 

and we have (see again [14]) 

dU E ,s 



(777.4.4) 



dUs 



ds 3 



(L'(0,L))' 



+ 



dUrs 



ds 3 



Ug\£ 



(L 2 (0,L)) 3 S 2 



ds 3 



US\£ 



<C 

(L 2 (0,L)) 3 5 



In order to be obtain the same estimate on Us that the one given by Lemma III. 3.1. in the previous section, 
we are lead to assume that \ug\s < CS K . Comparing with estimate (II. 3. 10) which gives 

\us\ £ <C(5 K + S 2K - 3 ) 

since ||dist(Va;«,5, SO(3))\\l2^ Vs - ) < CS K we are led to assume in the following that k > 3. 
The estimates (III.3.11), (III.4.1), (III.4.2) and (III.4.4) lead to the following lemma: 

Lemma III. 4.1. We assume that k > 3. There exists a subsequence (still indexed by S) of the sequence 
given in Lemma III. 3.1 such that 

1 3 
j^U^s — >U strongly in (77(0,7)) 

1 3 
j—[U E ,s^U E weakly in (i7 1 (0, L)) 



(777.4.5) 



^—^IZg^lZ weakly in (i7 1 (0,i))" 
(dUs 



Tight)- n Q 



\ds 3 

—ligus^-u weakly in (7 2 (0, L; H 1 ^))) 



Z a weakly in L 2 (0,L) 

3 



l^-T n 4V^ + (V^) T )-E' 
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weakly in (L 2 (f2)) 



3x3 



with E = (ni | n 2 | t)E (ni | n 2 | t) T where the symmetric matrix E is given by 



E 



/ du 1 r du du 

ds[- ni 2\ds;- n2 + ds- 2 

du 



dSi 



■ n 2 



V 



1 If du „ clTZ „ i \ 

1 r 95 „ dlZ _ i 
dft . „ dft . dt/ B 



-Si— ■ n 2 + S 2 — • ni 

ds 3 ds 3 



ds 3 



■ t 



Moreover the symmetric matrices E and E given in (III. 3. 10) satisfy 

1 



(777.4.6) 



E 




jli 3 -n. T n 



) if * 



if k > 3. 



Proof. The convergences (III. 4. 5) and the expression of E are proved in [14] taking into account (III. 4. 4) 
and the fact that \us\e < C5 K with k > 3. Let us notice that the limit U in the first convergence (III. 4. 5) is 
the same that in (III. 3. 3). In order to prove (III. 4. 6) we first write the Green-St Venant deformation tensor 
as (using u s = v s - Id) 



^tM(V x vs) T V x v 5 - I 3 ) = ((V x u s f + V x u s ) + i^Tls^iVxUsf^VxUs). 



Using convergences (III. 3. 5) and (III. 3.9), we deduce that 

, E' + -A T A if k = 3 
E= ( 2 

E' if s> 3, 

where the matrix A is defined in (III. 3. 3) and (III. 3. 5). Recalling relation (III. 3. 6) between A and TZ leads 
to (III.4.6). □ 

In the case k > 3, Lemma III. 4.1 shows first that starting from nonlinear deformations leads exactly 

to the same deformation model that starting from linearized deformations. The comparison in the case 

where k — 3 is more intricate. This is due to the definitions of the two warping v and u which do not 

satisfy the same geometrical conditions (see (II. 2. 7) for v and the beginning of this section for u). The 

second difference concerns the comparison between the stretching deformation Vs and the extentionnal 

displacement Ue- While it is easy to see that . t = ■ t for n > 3, in the case where n = 3 

«s 3 ds 3 



dU E dV s s 1 



dU 



one obtains — — • t = — — • t 

as 3 as 3 2 ds,^ 

contribution of the term tt^(R(5 — l3)t • (R$ — l3)t 



2 . 1 



. The correcting term 
2 6 2 



dU 



ds 3 



2 



actually comes from the limit 



IV. Asymptotic behavior of an elastic curved rod 

This section is devoted to use the above geometrical results in order to analyze the asymptotic behavior 
of an elastic rod when its thickness tends to 0. As usual, we consider a elastic energy density which is 
bounded below by dist 2 (F, 50(3)) (see e.g. [7], [11], [12], [19] and [20]). As mentioned in the introduction, 
our decomposition of the deformation permit us to scale the applied forces in order to obtain estimates on 
the deformation and on the total elastic energy (see (IV. 1.8)). Then, to simplify the argument, we specify 
the energy density through choosing a St Vcnant-Kirchhoff's material (see (IV. 1.9)). In the following, we 
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derive the limit elastic energy in the two cases k = 2 and k > 2 using T-convergence techniques. The limit 
energy is expressed as a functional of the fields V, R, V,s and v. Such limit energies depending on more 
variables than the 3D ones have also be derived in [19] by different techniques. In the present paper we also 
eliminate the fields Vs and v to be in a position to obtain a minimization problem for the rotation field 
R and the field V. Let us emphazise that in the r-limit procedure, the decomposition of the deformations 
is again helpful in two directions. Firstly it provides an explicit expression of the limit Green-St Venant 
deformation tensor and secondly it simplifies the proof of the two conditions involved in the identification of 
the limit energy by T- convergence. 

IV. 1 Assumption on the forces 

In this part we assume that the curved rod Vs is made of an elastic material. As in [5] and [11] we 
assume that the elastic energy W satisfy (actually we will consider an explicit energy) 

(JV.1.1) VF e M 3 , W(F) > Cdist 2 (F,50(3)), 

where C is a strictly positive constant. 

Let us denote by fs G (L 2 (Qs)) 3 the applied forces and by J{<j>) the total energy 

(7V.1.2) J{<j>) = f W{V4>) - f f s -4> 

This energy is considered over the set of admissible deformations: 

(JV.1.3) U s = {<j>e(H 1 (V s )) 3 \cj> = I d on T ,s}. 

For different boundary conditions see Subsection II. 2. 4. As far as the behavior of the forces fs is concerned 
we split the forces into two parts. The first one does not depends on the variables (si,S2) and the second 
part has a resultant equal to 0. Due to estimates (II. 3. 4), (II. 3. 6) and (II. 3. 7) the admissible order of theses 
forces can be chosen different. 

Let / be in (L 2 (0,L)) 3 and let g be in (L 2 (Q,)) 3 such that 

(IV.IA) f g(S 1 ,S 2 ,s 3 )dS 1 dS 2 = for a.e. s 3 e]0,L[. 

J UJ 

We assume that fs is denned by 

(IV.1.5) f 5 ( s ) = 6 "f(s 3 )+5 K - 1 g(^-,?f,s 3 ) fora.e. seSl, 

The fact remains that to find a minimizer or to find a deformation that approaches the minimizer of J((j>) 
or of J(<t>) — J(Id) is the same. Let v be in U$, thanks to (II. 3. 7), (II. 3. 9), (IV. 1.4) and (IV.1.5), we obtain 

(IV.1.6) f h-{v-I d ) <C5 K (||/ (L2(0 . L))3 + || ff || (i2(o)) 3)||dist(V«,5O(3))|| i2(n) . 

Actually one can think to use estimate (II. 3. 4) instead of (II. 3. 7) and (II. 3. 9) in the above inequality. The 
reader will easily see that this gives a better estimate only in the case where n < 2 which is not considered 
in the following (see Remark below). 
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It is well known that generally a minimizer of J does not exist on XJs- In the next sections we will 
investigate the behavior of the functional -^{Jfy) — J {Id)) m the framework of the T-convergence. Hence, 
we assume that 

(IV.1.7) _L (./(„) _,/(/„)) <d 

where C\ does not depend on 8. Using (IV. 1.1) and (IV. 1.6) we obtain for such v 

C||dist(V^, S'OCS))!^,^^ - C^(||/|| (i2(0)) 3 + ||.9|| (i2(f , )) 3)||dist(V^, S'0(3))|U 2(7 , ts) < d5 2 -. 
Hence, we have 

(JV.1.8) ||dist(V«,SO(3))|| i2(7 , 4) < C5 K . 

where the constant C depends on the sum ||/||(L 2 (o,i)) 3 + ll</ll(L 2 (fi)) 3 and of C\. 

Remark. If one uses estimates (II. 3. 4) to bound the contribution of the forces in the energy, one alternatively 
obtains through similar calculations 

||dist(Vi;,50(3))|| L2( ^ ) <C ( 5 1 + K / 2 . 

Comparing to (IV. 1.8), one gets a better estimate only if k < 2. 

Let us notice that once the assumption (IV. 1.5) on the applied forces is adopted, the estimate (IV. 1.8) 
and the results of Section II permit to obtain estimates of V, R, v and V x v — R with respect to 5. To 
emphasize how these estimates can help pass to the limit as S tends to 0, we will restrict the following 
analysis to a classical and simple elastic energy. We denote by tr{A) the sum of the elements on the main 
diagonal of the 3x3 matrix A. 

In order to simplify the derivation of the limit model we choose 

(IV.1.9) W( F)=(^^ FTF - l3 )) 2 + ^^ TF - 13 ) 2 ) * det(F)> ° 

[ + oo if dct(F) < 

which corresponds to a St Venant-Kirchhoff 's material (see [7], [8]). 
For all 3 x 3 matrices such that det(F) > we have 

W{F) > J|||F T F-I 3 ||| 2 > ^dist 2 (F,50(3)) 

hence assumption (IV. 1.1) is satisfied. For every ^eUj satisfying (IV.1.7), we have using (IV. 1.6) 

^||V0 T v^ - h\\ 2 {LHPs))3x3 < J(4>) - J{h) + [ fs-{4>- id) 

< C 1( 5 2 « + C^^(||./-|lci 2C o,i))3 -f- ||^|| c ^ 2C ^ )) 3)|| c list(V0, 5'0(3))|U 2C ^ 5) . 
Due to estimate (IV. 1.8) we obtain the following estimate of the Green-St Venant's tensor: 

{IV.l.lQ) ll^{W T V0-I 3 }|| (L2(n)) 3 X 3 <C6". 
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It results that <f) belongs to (W 1,4 (Ps)) 3 and moreover 



(JV.1.11) 



I|V0||(L4(^))3X3 < C(5 5 . 



Furthermore, there exists two strictly positive constants c and C which does not depend on S such that for 
any <j> £ Us satisfying (IV. 1.7) we have 



(J V. 1.12) 



-c8 2K < J(<j)) - J{I d ) < CS 2K . 



We set 



m s = inf (J(<j>)-J(l d )). 



(JV.1.13) 

As a consequence of the inequality in (IV. 1.1 2) we have 
(JV.1.14) 
We denote 



(JV.1.15) 



m K = lim — . 

<5-vo S 2K 



IV. 2 Limit model in the case n = 2 

Let ( w <5) < < 5<(5 be a sequence of deformations belonging to U,5 and such that 



(JV.2.1) 



,. J(vs)-J(Id) , , 
lim < +oo. 

5^0 <5 4 



Upon extracting a subsequence (still indexed by 5) we can assume that the sequence (vs) satisfies the 
condition (IV. 1.7). From the estimates of the previous section we obtain 



(iV.2.3) 



f ||dist(V^,50(3))|| L2(n) <CS 2 , 
lli^vJVvs -I 3 }\\ {L2{Vs))3x3 

k ||V« 4 || (L 4 (7 , 4)) 8 X S < CS^ 



<CS 2 



For any fixed S € (0, Jo], the deformation vs is decomposed following (II. 2.1) in such a way that Theorem 
II. 2. 2 is satisfied. There exists a subsequence still indexed by 6 such that (see Section II. 6) 



(JV.2.4) 



3x3 



' Us R weakly in (tf^O, L)) 
V 4 — -> V strongly in (if^O, L)) 3 
Vb.S^V strongly in (H\<d,L)f 

\v s ,s^Vs weakly in {H\U,L)f 
^U 5 v 5 ^v weakly in (L 2 (0, L; H l (w))f 
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where R(s 3 ) belongs to 50(3) for any s 3 <= [0,L], V e (# 2 (0,L)) 3 together with 



(JV.2.5) V(0) - M(0), R(0) = I 3 , V s (0) = 0, and — = Rt. 

ds 3 

Furthermore, we also have (see (III. 2. 6) and (III. 2. 7)) 

'Tlsvs—tV strongly in (^(fi)) 3 , 



(JV.2.6) 



where 



v 3x3 



n 5 (V xVl5 ) — ► R strongly in (L 2 (fi)) ; 
-^n 5 ((V x ^) T V x ^-I 3 ) weakly in (L 2 (tt)) 3x3 , 



E 



-i{NI-l*>(Sl 



— |-(5 ini + 5 2 n 2 ) + — ) R 



Remark that, due to the decomposition (II. 2.1), the convergences (IV. 2. 4) and (IV. 2. 6) imply that 
(IV.2.7) 



U5{V5 X Vs) — >Si(R- 13)11: + S 2 (R-I 3 )n 2 strongly in (L 2 (Ct)) 3 . 



Now, recall that 

( J(v 5 ) - J{I d ) 



<5 4 



(JV.2.8) < 



{ S K ^ n *(( v ^«) Tv ^« - I a))] 2 + iA 1 1 ^n 5 ((v^ 5 ) T v^ 5 - 1 3 )|| | 2 }|n 5 dct(v$)| 



- [ f-iL s (v s -i d )\n 5 det(vn- I g Us{v * /d) |n^det(v$)| 

v Jn Jn 



In order to obtain the limit of the terms involving the forces, we recall that det(V$) = 1+Si det ( rii I n 2 I — + 

V as 3 / 

s 2 dct ^ni I n 2 | ^r~~) so that indeed H5 det(V$) strongly converges to 1 in L°°(Q) as 5 tends to 0. As a 
consequence and using the convergences (IV. 2. 6) and (IV.2.7), it follows that 

Rsivs - Id) 1 



8^>Q 
r-L 



lim / / 



n 4 (««-/ d )|IIadet(V$)|+ I 9 

Jn 



\Us det(V*)|) 



= £ (\"\f + 9S a det (m|n a I ^)d5idS 2 ) ■ (V - M) + ^2^ (J gS a dS 1 dS 2 ) -(R-I 3 )n Q . 

In order to pass to the limit-inf in the left hand side of (IV. 2. 8), we recall that the map M 1 — > ^(tr(M)) 2 + 

y(i|||M||| 2 is continuous and convex from M 3 into M, so that the map A 1 — > I (^(ir(^4)) 2 + ,u||m|| 2 ) from 

(L 2 (f2)) 3x3 into R is weak lower semi-continuous. The above strong convergence of IT,5det(V$) together 
with convergences (IV. 2. 6) finally give 

' /^(ME)) 2 + ,|||E||| 2 } 
(/V.2.9) < ~ £ (Hf + J pJjS a det{n 1 \n 2 \^)dS 1 dS 2 )-(V-M) 

- V / { I gSadStdS^) ■ (R - I 3 )n Q < lim ^(J(v s ) - J{h))- 
~t J a V Jlj ' <5^0 
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Let XJniin be the set 

U nKn = {(V',R',V;,TJ') e (H 2 (0,L)f x (^(0,L)) 3x3 x (H\0,L)f x (L 2 (0, L; H 1 ^))) 3 \ 

V'(0) = M(0), Vg(0)=0, R'(0)=I 3 , f v'{S 1 ,S 2 ,s 3 )dS 1 dS 2 = for a.e. s 3 e (0, L) 

R'(s 3 ) € 50(3) for any s 3 € [0,L], — = R't}. 

as 3 J 

The set U n iin is closed in the product space. For any (V ,R ,V s ,v) G U n lin, we denote by Jnl the 
following limit energy 

' J NL {v',R,V' s ,v') = jf {^(tr(E)) 2 + mIIIE'IH 2 } 
(JV.2.10) < - 1 (nf + ^JgS a det(n 1 \n 2 \^)dS 1 dS 2 )-(v'-M) 

- E l L ( / gSadStdSi) ■ (R - I 3 )n Q 



where 



(JV.2.11) 



. +R 1 a& 1 d^ 51nl + 52n2) + *? ) (nl 1 n2 1 1} }■ 

With this notation (IV. 2. 9) reads as 

(/V.2.12) ^(V,R,V S ,U) < Um ~ -TO 

,5^0 <5 

Now let (V,R,Vs,w) be in U^n and let ^(Vjv, Rjv, Vs,n,vn)^J be a sequence of elements belonging 
to Unlin such that 

V N e {W 2 '°°{0,L)f, V N ^V strongly in (H 2 (0,L)f 
Rjv e (W rl ' oo (0,L)) 3x3 , Rw — >R strongly in (JT^O, L)) 3x3 
(/r.2.1:-!) <( Vs.atG (W 1 ' oo (0,L)) 3 , Vsjv^Vs strongly in (ff^L)) 3 

Fat g (W 1>0 °(ft)) 3 , v N (S u S 2 ,0) = 0, fora.e. (5i,5 2 )gcj, 

Uat — >v strongly in (L 2 (0, L; H 1 ^)) 3 . 

To prove the existence of the sequence (Rjv) .wen* see the appendix at the end of the paper. 
We consider the deformations (S € (0,<5o]) 

{IV.2.U) v N , s (s) = V N (s 3 ) +K N {s 3 )(s 1 n 1 + s 2 n 2 ) +SV s ,n{s 3 ) +S 2 v N ^,j,s 3 y sett s . 

Using convergences (IV. 2. 13), the fact that (Vjv, Rjv, V.s,n,vn) belongs to U nJi „ and proceeding as in Sub- 
section III. 2, we have 



(JV.2.15) 



HsVn,s—>Vn strongly in (W 1 ' 00 ^)) 3 , 
U s (V x v NtS ) — > Rjv strongly in ( J L°°(f7)) 3x3 , 



—U 5 ((\7 x v N ^) T ^ x v N ^ - I 3 ) 



E 



N 



strongly in (L°° (fl)) 3x3 , 
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where 

i7 If/ I , .fdv N dv N dR N , rfV s ,w\ T „ 
Ejv =-{ (ni | n 2 1 | — | -^-(S 1 n 1 + S 2 n 2 ) + -^-J 

If 5 is sufficiently small we have det (V x vn,s(x)) > for a.e. x £ Vs because of the second convergence in 
(IV. 2. 15). It follows that J(vjy^) < +oo. We divide J(vn,s) — J(Id) by 8 4 and we pass to the limit using 
the strong convergences (IV. 2. 15). We obtain 

(iV.2.16) lim-U^K,*)- J(I d )) = Jnl(Vn,Rn,V n . s ,vn)- 

Letting TV tend to +oo and using (IV. 2. 13), it follows that for any (V, R, Vs, v) € U n / in 
(/V.2.17) Jnl(V,R,V s ,v) = N ma jo J NL (VN,B. N ,V N ,s,VN). 

Hence, through a standard diagonal process for any (V, R, Vs, v) € U n lin there exists a sequence of admissible 
deformations vs & (H 1 (Vs)) 3 such that 

(/V.2.18) Jnl(V,R,Vs,v) = lim JM ~ J[Id) . 

(5^0 

The following theorem summarizes the above results. 

Theorem IV. 2.1. The functional J NL is the T-limit of "^'^ - ^ ^ e following sense: 

• consider any sequence of deformations ( v s) 0<s<Sg belonging to Us and satisfying 

y J(v S )-J(Id) „ 

lim -7 < +oo 

and let (Vs, R<5, Vs t a, 1^5) &e f/ie £erms 0/ i/ie decomposition of vs given by Theorem II. 2. 2. Then there exists 
(V,R, Vs,v) G Unlin such that (up to a subsequence ) 

Rj^R weakly in (H 1 ^, L)) 3x3 
Vs — >V strongly in (H^O^L)) 3 
Vb.S — > V strongly in (iJ^O,^)) 3 
^V S , S ^V S weakly in (^(O.L)) 3 

II5W5 — w weakly in (L 2 (0, L; i? 1 (w))) 



<5 2 



and we have 



J NL (V,R,Vs,v)<M J{Vt) - J{Id) 



5^0 

/or any (V,R, Vs,u) G U„;,n i^ere errasis a sequence ( v s) 0<s<s belonging to XJssuch that 



Jnl(V,R,V s ,v) = lim 



0<5<5 
J(^) - J(Id) 



8^0 <5 4 
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Moreover, there exists (Vo, Ro, Vs,o, vo) € U„; in sucft i/iai 



(JV.2.19) 



m 2 = lim — = Jjvl(V ,Ro, Vs,o, «o) = mm Jnl(V, R, Vs,tJ). 

" (v,R,Vs,w)eu n!in 



The next theorem shows that the variables Vs and v can be eliminated in the minimization problem 
(IV. 2. 19). To this end let us first introduce a few notations. We denote by E the Young's modulus of the 
material and by \ the solution of the following torsion problem: 



(/y.2.22) 



At least we set 



k = 



95 2 



•Si 









r: 


/i = / si 


'■<=! 




J UJ 


J UJ 



si 



Theorem IV. 2. 2 Let (Vo,Ro) be given by Theorem IV. 2.1. The minimum mi of the functional Jnl over 
U„Hn satisfies the following minimization problem: 



(JV.2.20) 



where 



m 2 = Fnl(Vq,Hq) = min F NL (V,H), 
(v,r)gv„ h „ 



V nlm = {(V,R) G (H 2 (0,L)) 3 x (^(O.L)) 3 * 3 | V(0) = M(0), R(0) = I 3 , 
R(s 3 ) G SO(3) for anys 3 G [0,L], ^ = Rt), 



and 



(JV.2.21) < 



- / L (l w l/ + E / 5^det(m |n 2 | ^)dSidS 2 ) ■ (V - M) 
-Y. t { I 9SadS 1 dS 2 ) ■ (R - I 3 )n a . 



Proof of Theorem IV.2.2. Let us first notice that in the expression (IV.2.10) of J NL (V,TL,Vs,v), one 
can replace E by E where E and E are given by (III.2.7), (III.2.8) and (III.2.9). 

In order to eliminate (Vs,v), wc fix (V, R) G V n lin and we minimize the functional Jnl(V, R, ■, •) over 
the space 

W={(V s ,v)e(H 1 (0,L)) 3 x(L*(0,L-,H\u>))) 3 | Vs(0) = 0, 
f v{S 1 ,S 2 ,s 3 )dS 1 dS 2 = for a.e. s 3 G (0, L)\. 
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Through solving simple variational problems (see [14] ) , we find that the minimum of the functional Jn l (V, R, • , •) 

over the space W is obtained with ■ Rt = and 

ds 3 



(/V.2.23) 



2 ds^ 



t • Rm + 5i&— t-Rn4 
ds 3 ) 

S2 — Si dR 



v{S 1 ,S 2 ,.)-B.n 1 =i/{ 

v(S u S 2 , •) • Rn 2 = vfs&^t ' Rn i 
I- ds 3 

v(Si,S2,0-ra + Si^-Rni+S 2 ^-Rn 2 = {xO^, S 2 )^ru • Rn 2 }Rt 

as 3 as 3 ds 3 



t -Rn 2 } 



where 1/ 



A 



is the Poisson's coefficient of the material. Then the symmetric tensor E (see again 

2(A + /i) 

(III. 2. 9)) at the minimum is given by 



/ -^E 33 (R) 



(7V.2.24) 



E(R) 



~ . 1 / dy „ \ dJU 
-i/E 33 (R) -(^ 

* E 33 (R) 



ni • Rn 2 



/ 



rfR 



dR 



where E 33 (R) = —Si- — t • Rni — S2- — t • Rn 2 . Upon replacing E by E(R) in the expression of Jnl we 



obtain 



ds 3 ds 3 

min J NL (V, R, V s ,v)= Fnl (V, R) , 
(v s ,«)ew 

where the functional Fnl is given by (IV. 2. 21). 

Remark. The above analysis shows that if ( v s) 0<s<So is a sequence such that 

J(vs) - J(I d ) 



□ 



m 2 = lim 

<5->o 



<5 4 



then there exists a subsequence and (Vo,Ro) € V n u n , which is a solution of Problem (IV. 2. 20), such that 
the sequence of the Green-St Venant's deformation tensors satisfies 



^U 5 ((V x v 5 ) T V x v 5 - I3) — > (m I n 2 I t)E(Ro)(ni | n 2 | t) T strongly in (L 2 (0)) 3x3 , 

where E(R ) is defined in (IV.2.24). 
IV.3 Limit model in the case k > 2 

be a sequence of deformations belonging to XJg and such that 



Let H <S<So 

(IV.3.1) 



lim — < +00. 



.5^0 



Upon extracting a subsequence (still indexed by 5) we can assume that the sequence (v$) satisfies the 
condition (IV. 1.7). From the estimates of the section IV. 1 we obtain 



{IV.3.3) 



( \\dist(\7vs,SO(3))\\ L 2 (P5) <CS K , 

I ll^{v^v^-i 3 }|| (i2(n))3>! 3 <CS K , 

[ l|Vu5||(i4 (n)) 3x3 < CSi. 
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For any fixed S G (0, So], the displacement us = v$ — Id is decomposed following (II. 2.1) and (III. 3. 2) in such 
a way that Theorem II. 2. 2 is satisfied. There exists a subsequence still indexed by S such that (see Section 
III.3) 



(JV.3.4) 



^(R«-l3)^A weakly in {H\0,L)) 
-L^Us^U strongly in (H\0,L)) 3 



3x3 



-Vs,8^V s weakly in (H\0, L)f 



-Ugvs^v weakly in (L 2 (0, L; H 1 (lj))Y 



where A is an antisymetric matrix and U(Qi) = Vs(0) = 0. Moreover, U belongs to (H 2 (Q,L)) 3 and there 
exists Tl G (H 1 ^, L)) 3 with 11(0) = such that 



(JV.3.5) 



Furthermore, we also have 



^- = At = K A t. 

ds 3 



r 1 



S 



(JV.3.6) 



^n 5 u,5 — >W strongly in (if 1 (ft))' 



n 5 (v x ^) — 
n^(Mg - Z4) 



I 3 strongly in (L 2 (fl)) 3x3 , 
— > S-JL A ni + #2 ft A n 2 strongly in (i 2 (ft)) 3 , 



^ — T U s ((V x vs) T V x v s — I3) — E weakly in (L 2 (ft)) 3x3 , 
where (see Section III. 3) 

E = (ni I n 2 |t)E(ni | n 2 |t) T 



(JV.3.7) 



/ dv 1 ( dv dv 1 If 9» _ dTZ dVs 1 \ 



E 



^2" 



n 2 



2W t + s W t+ d^' n2 } 

-Si- n 2 + S^-p- • ni + — — • t , 

ds 3 ds 3 ds 3 / 



Proceeding as in the previous section, we pass to the limit-inf in IM^M and 

we obtain 



S 2k 



(IV.3.8) 



+ HIIEIH 2 } 

2 



f/„{> <E » 2 

- I 9(Si,S 2 ,.)S a det(m I n 2 I ^d&dSa) U 

- V [ L ( f g(S 1 ,S 2 ,.)S a dSidS 2 ) ■(RAn„) < Vm-L(j(v s )- J(I d )). 
a= iJo y Ju> ' W 



Let Ujin be the space 



U K „ = {{U' ,Tl' ,V' s ,v') G (# 2 (0,L)) 3 x (H^O^jf x (H\0,L)) 3 x (L 2 (0, L; H 1 ^))) 3 



z/(0) =ft'(0) = Vo(0) = 0, ^-=KAt, 

ds 3 



{ v'{S 1 ,S 2 ,s 3 )dSidS 2 = for a.e. s 3 G (0,L)}. 
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For any (ll , 1Z , V s , v ) e U; in , we set 

' j L (u',n',v' s ,v') = J a {^(tr(E')) 2 + hhe'ih 2 } 

- £ (Wf + YlJ g(Si,S 2 ,.)S a det (m |n 2 | ^)dS 1 dS 2 ) -U 



(IV.3.9) 



f L ( / <nS t ..S. 2 . ;.s'.. ,/>':</>,) ■ (7? u„). 

», a=l ' 



with E is given by (IV. 3. 7) where we have replaced (U, 1Z,Vs,v) by (U ,TZ ,V s ,v ) . From (IV. 3. 8) it results 
that 



(IV.3.10) 



J L (U,1l,Vs,v) < Ijm ipi 



.^J(v s )-J(I d ) 

5^0 



Now, let (U, 1Z, Vs,v) be in TJu n and let [IZn, Vs n,vn ) sequences of elements such that 



(JV.3.11) 



Moreover we set 



K N (0) = 0, Vs,j V (0) = ) v N (S 1 ,S 2 ,0) = 0, for a.e. (S 1 ,S 2 )eu, 
TZn € (T4^ 1,oo (0, L)) 3 , TZ N ^TZ strongly in (H 1 ^,^) 3 
V s , W e (VF 1 ' 00 ^)) 3 , V s . N ^Vs strongly in (ff^L)) 3 



v N G (W 1,00 (fi)) 3 , U;v — >TJ strongly in (L 2 (0, L; H^uj)) 



iV 



ds 3 



= ^ArAt, Wjv(0) = 0). 



We consider the deformations (S € (0, Jo] and s € ^,5) 
(/V.3.12) v N , 5 (s) = V N ,s(s 3 ) + Riv, 4 (*3)(*ini + s 2 n 2 ) + S^VsM^) + ^(j- J> s z) 

where Rn,s and Vn,s are defined below 



dR NiS _ c K -2 R id 



d,S3 

RiV,5(0) = i 3 



/■S3 

Vat. 5 (s 3 ) = M(0) + / R NtS (z)t(z)dz. 
Jo 

Here B is the 3x3 antisymmetric matrix such that 

VxeM 3 , B N x = — — A x. 

ds 3 

Using the above convergences and the fact that (Un,TZn,Vs,n,vn) belongs to Uu n , we have 



(JV.3.13) 



— ^U s u NtS ^U N strongly in (VF 1 ' 00 ^)) 3 , 
n 5 (V x Wjv,5) — > I 3 strongly in (L°°(Q)) 3><3 . 
^ ^^((V^afV^iv.a - I 3 ) — ► strongly in (L°°(tt)) 3x3 , 
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where E^v is given by (IV. 3. 7) where we have replaced (l{,TZ,Vs,v) by (Un,1Zn, Vs,n,vn)- If S is sufficiently 
small we have det (V x vn,s{x)) > for a.c. x G Vs- We divide J(vN t s) — J{Id) by S 2k and we pass to the 
limit. We obtain 

(/V.3.14) lim -L(J(^ >5 ) - J(/ d )) = J l (Un,Kn,Vs,n,vn). 

Now letting N tend to +oo gives that for any (U,1Z, Vs,v) G Uu n 

(/V.3.15) Jl(U,K,Vs,v) = lim Jl(U n ,K n ,V s ,n,v n ). 

JV->+CX> 

Hence, through a standard diagonal process for any (U, 1Z, Vs,v) G Ujj n there exists a sequence of admissible 
deformations vg G (H 1 (Vs)) 3 such that 

(/V.3.16) J £ (M,ft,y s ,P) = lim IMJZW. . 

The following theorem summarizes the results of the case k > 2. 

Theorem IV. 3.1. T/ie functional Jl is the T -limit of "^'^ ^l^^ in the following sense: 
• for any sequence of deformations ( v s) 0<s<Sg belonging to Us and satisfying 

hm < +oo 

and let (lAs,TLs, Vs,6,vs) be the terms of the decomposition of the displacement us = vs—Id given by (III. 3. 2). 
Up to a subsequence there exists (l4,lZ,Vs,v) G Uu n such that 



-^-L(R 5 -I 3 )^A weakly in (H\0, L)) 3x3 
J^_ 2 Us — >U strongly in (_ff 1 (0,L)) 3 
J^Vs.s^Vs weakly in (^(O.L)) 3 

1 3 

—Hsvs^v weakly in (L 2 (0, L; 



(5 

where for any x G M 3 , Ax — TZ A x and we have 

j L (u,n,v s ,v)<m J{v5) - J{Id) 

s^ S 2K 

• for any (U,lZ,Vs,v) G U^n there exists a sequence ( v s) 0< s<s belonging to XJ$ such that 

J L (U,ll,Vs,v)=^ J ^- J ^\ 

Moreover, there exists (l4o,7lo 7 Vs,o,vo) G Uu n such that 

m « = ] im = ^l(Uo,T^o,Vs,o,vo) = min J L (U,TZ,V s ,v). 
° {u,n,Vs,v)ev lin 
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The next theorem is the analog of Theorem IV. 2. 2. 

Theorem IV. 3. 2 Let (Uo,7Zo) be given by Theorem IV. 3.1. The minimum m K of the functional Jl over 
Uhm satisfies the following minimization problem which admits a unique solution: 



(JV.3.17) 



where 



and 



m K = T L (U , 7Z ) = min T L (U, U) 



V Kn = {(U,7l) e (H 2 (0,L)f x (ff^O.L)) 3 | H(0) = ft(0) = 0, £-7iAt} 



(7V.3.18) 



r L (u, K) = -l (_-„ 2 ) + _y o ( 



ds 3 



/lis: 



/dTZ \2 



(W + E / g(Si,S 2 ,.)S a det(n 1 | n 2 | J^)dS 1 dS 2 ) . U 
-V [ L ( [ 9 (S 1 ,S 2 ,.)S a dS 1 dS 2 ) .(RAnJ 



-E is i/ie Young's modulus and K is given in Theorem IV. 1.3. 

Proof of Theorem IV.3.2. We proceed as in Theorem IV. 2. 2. We fix (V, 1Z) e Vu n and we minimize the 
functional Jl(U,TZ, ■, •) over the space W. Through solving simple variational problems (see [14] again), we 

find that the minimum of the functional jAlA,TZ, ■, ■) over the space W is obtained with • t = and 

v ' ds 3 



(JV.3.19) 



v(Si, S 2 , ■) • ni = -z/j 



Sf - 5? dTZ „ „ dK ~\ 

n 2 + S 1 S 2 —-n 1 \ 

ds 3 j 



2 ds 3 

,n r, ^ { r, r, dll S% ~ S? dTZ 

v{S u S 2 ,.)-n 2 = -v{-S 1 S 2 —-n 2 + ~ 



2 ds 3 



-} 



S 2 , .) • t + • m + • n 2 = X (S U S 2 )^ ■ t 

as 3 as 3 ds 3 



Then the symmetric tensor E (see (IV. 3. 7)) at the minimum is given by 

* * E 33 (7e) 



(JV.3.20) 



E(K) = 



V 



where E 33 (7£) = — 5i4— ■ n 2 + S2-7— ■ ni. Upon replacing E by E(7£) in the expression of J7l we obtain 

ds 3 cfs 3 

min J L (U, 71, V s ,v)= T L (U, 71) 
(v s ,w)ew 

where the functional JFl is given by (IV. 3. 18). 

Remark. The above analysis shows that if ( v s) 0<s<So is a sequence such that 

J(vs) - J(Id) 



□ 



m K = lim 



■5^0 8 2k 
32 



then the sequence of the Green-St Venant's deformation tensors satisfies 

y^M^*v s ) T V x v s - I 3 ) — ► (m I n 2 I t)E(72o)(ni | n 2 | t) T strongly in (L 2 (0)) 3x3 , 

where E(ft ) is defined in (IV.3.20) and K is the solution of (IV.3.17). 
IV.4 Extentional models for special forces. 

In this subsection we investigate the case where f$ is given by 

(/V.4.1) f s (s) = S*-\f(s 3 ) for a.c. a E fi 4 , 

where / belongs to (L 2 (0, L)) 3 . Without any additional assumption on /, Subsection IV. 1 shows that this 
leads to 

||dist(V« 5 ,SO(3))|| L2(n) <CS K -\ 

if ^ 2 K-2 (J( Vs ) ~ J(Id)) < Ci- As a consequence, the results of Subsections IV. 2 et IV. 3 can be applied if 
k > 3). Let us for example consider the case k > 3 and remark that due to the choice of / the contribution 
of the forces in the limit energy Fl{Uq,TZq) is equal to 

-MjT f{s 3 )-Uo{s 3 )ds 3 = -\u\ /(s)ds) -^o(s 3 ) At(s 3 )ds 3 . 

Then if the quantity / f(s)ds is proportionnal to t(s 3 ), this contribution vanishes and then IZo = Uq = 

J S3 

and the minimum is null. This example shows that for this kind of special forces, the energy have a smaller 
order than 2k — 2 or equivalently that the estimates on vg can be improved in this case. 

We assume that there exists / e _ff 1 (0,i) such that 

(IVA.2) f f(l)dl = f(s 3 )t(s 3 ) for any s 3 e [0,L\. 

J S3 

Let v an admissible deformation of the rod Vs- Now, using (IV. 4. 2) we derive a new estimate of 

/ /«■(«- Id)- Notice first that det(V$) = 1 + s 1 dct (rii | n 2 | — ) + s 2 dct (ni | n 2 | — V then 
jp { V ds 3 ) \ ds 3 ) 

using the decomposition (II. 2.1) for the admissible deformation v, estimates of Theorem II. 2. 2 and (II. 3. 5) 
together with / s a ds\ds2 — we deduce that 



(JV.4.3) 



/ fs-{v-I d )-\uj\5 K+l [ f(s 3 )-(V(s 3 )-M(s 3 ))ds 3 

JVs JO 

l<^ +1 ||.f|| (i2( o,L))3||dist(V«,50(3))|| i2(n) . 



We obtain by integrating by parts and using the decomposition (II. 2. 16) of V (see also (III. 4. 7)) 



(jy.4.4) 



J /(*») ■ (V(« 3 ) - M( S3 ))ds 3 = J /(« 3 )t(*3)-(^(*3)-t(s 3 ))d83 

: / L /(*3)t(* 3 ) ' (R(«3) - I 3 )t(s 3 )dS3 + / L /(« 3 )t(« 3 ) ' ^(S3)^3 
JO JO « s 3 



ds 3 
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= - J /( S 3)(R(S 3 ) - l3)t(* 3 ) • (R(«3) - l3)t(* 3 )d*3 + /(*3)^(*3 ) • t(.s :i lr/.s,. 

I ^ Jo Jo 



Finally, using the above estimate and (II. 2. 17) we get 



(JV.4.5) 



We assume that 



/ fs ■ (V - Id) + / /(«3)(R(« 3 ) - Is)t(s 3 ) • (R(«3) - Is)t(« 3 )d«3 

JVs 1 JO 



V 5 

l<^ K ||/|| (i2( o, L) )3||dist(V U ,50(3))|| L2( ^ ) . 



IVs JVs 

which implies using (IV. 1.9) 



/ W(Vv) - [ f 5 -(v- I d ) = J(v) - J(I d ) < +^ 

JVs JVs 



(/V.4.6) £||dist(Vi;,SO(3))||£ a(7 , 4) - / fs-(v-I d )< J(v) - J(I d ) < +oo. 

4 JVs 



Hence 



(IVA.7) 



£||dist(Vt>, SO(3))\\l 2{Vs) + J /(R - I 3 )t • (R - I 3 )t 

< CS K \\f\\ {L 2 {QiL)yi \\di S t(\7v,SO(3))\\ L 2 {P5) + J(v) - J(I d ) < +^ 



Now, in view of the above inequality, let us consider a sequence ( v s) 0<s<So satisfiyng J(v$) — J(I d ) < C\5 2k . 
From estimate (II. 3. 8) we deduce that 

(7V.4.8) /^/(R-IsJt^R-IsJt <C*<5- 4 ||/|| (i2(0 . L)) 3||dist(V U ,5O(3))||i 2(n) , 

Jo 

where the constant C* only depends upon the geometry of the middle line of the rod. According to inequalities 
(IV. 4. 8) and (IV. 4. 7) the cases n = 3 and k > 3 lead two different energy estimates. 
If k > 3 we obtain 



(JV.4.9) 



' ||dist(V^,50(3))|U 2(n) <CS K , 
\\^{WjVv s - I 3 }|| (i2(n))3X 3 < CS K , 



The constant C does not depend on S. 

If k — 3, the energy estimate depends on 1 1/| | (x, 2 (o.z,)) 3 - Indeed , if 
(IVA.10) 11/11(^(0,^)3 < 



2C*\w\ 

estimate (IV.4.8) and (IV.4.7) give 



M iiji|(i,2(o 



, i)) 3)||dist(V^,50(3))||| 2(n 



(7V.4.11) { V4 2 

< C<5 3 ||/|| (L2(0 , L)) 3||dist(V«,5O(3))|| L2(n) +C7 1( 5 6 



and then 



f ||dist(V U5 ,50(3))|| i2(n) <C5 3 



(JV.4.12) <j ||i{V W JV^-I 3 }|| (L2( ^ ))3X 3 <C5 3 , 



l|V« 4 || (L 4 (7J4)) SXS < C(5J 
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The constant C does not depend on 5. 

In view of (IV. 4. 7), an alternative assumption to (IV. 4. 10) which also leads to estimates (IV. 4. 12) is to 
suppose that 

(7V.4.13) f(s 3 ) > for almost any s 3 e{0,L). 

In the sequel of this subsection, we will assume that the forces f satisfy (IV. 4-2) together with (IV. 4-10) 
or (IV. 4. 13) ifn = 3. 

Now we have to pass to the limit-inf in — ^ g 2l f^ ^ ■ According to estimates (IV. 4. 9) and (IV. 4. 12), 
performing this process in the elastic energy term is identical to the one detailed in the previous Subsection. 
We just focus on the behavior of the terms involving the forces. In view of (IV. 4. 3), (IV. 4. 9) and (IV. 4. 12) 
we get 

iTo W J Vs h ' (VS h) = JTo 6^ £ /(S3) ' (Vi(53) _ M ^ ds3 - 
Now we use the notations and results of Subsection III. 4, we have 



Thanks to the convergences of Lemma III. 4.1 we deduce that 
(IVA.U) }z ^j^f 5 . { v 6 -I d ) = Hjjt. 

Let us define the limit operator Jls by 



cIUe.s 
ds 3 



dU E 
ds 3 ' 



(IVAA5) V(U,1l,Vs,v)eUi in , Jls(U,K,Vs,v) - ^ {^(tr(E)) 2 + M |||E||| 2 } - | w | ^ / ^ • t. 

The matrix E is given by (IV. 3. 7) and the displacement Ue is such that (see Lemma III. 4.1) 

dU 



dVs t 1 
dUs A J ds 3 2 



" , ' J -"" si-' i .a 

ds 3 



ds 3 

• t if k > 3. 



if k = 3, 



The expression of Jls shows that this functional has a unique minimizer. We have obtained the following 
result. 

Theorem IV.4.1. The functional J LS is the T-limit of pt^^ * n ^ e f°H owin 9 sense: 
• for any sequence of deformations ( v s) 0<s<So belonging to TJ$ and satisfying 

,. J(vs)-J(Id) „ 
km 7 < +oo 

,5^0 S lK 
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and let (lis, R<5, Vs.s, v s) be the terms of the decomposition of the displacement us = vs — Id given by (III. 3. 2) 
and (III. 4-3). Up to a subsequence there exists (jJ,TZ,Vs,v) € \5u n such that 



.3x3 







1 ' 




5 K - 




1 




S K ~ 




I— 1 





1 3 

— Il^tJ^ — U weakly in (_L 2 (0, i; 1 (w))) 

w/iere /or any x e R 3 , As — TZAx and where the relation between U E , Vs andU is given by (IV. 4.16). We 
have 

J(v s ) - J(I d ) 



J LS (U,n,V s ,v) < lim 



<5-s-0 



6 2k 



• for any (U,7Z,Vs,v) € U ;i „ there exists a sequence ( v s) 0<s<So belonging to XJssuch that 

Moreover, there exists a unique (Uo,Tto, Vs,o,^o) € U/ in such that 

m « = ] im = <?ls(Uo,Ko,Vs,o,vo) = inf J LS (U,K,V s ,v). 

(u,n,Vs,v)ev Hn 

The next theorem is the analog of Theorems IV. 2. 2 and IV. 3. 2. 

Theorem IV.4.2 Let (U ,K ,Vs,o) be given by Theorem IV. 4.1 and U E ,o € D Ex defined by (IV. 4.16). The 
minimum m K of the functional Jls over U; in is obtained with Uq = TZq = and it is given by the following 
minimization problem which admits a unique solution: 

(/V.4.17) m K = F LS (U E fi) = min F LS {U E ), 

where 

E is the Young's modulus. 

Proof of Theorem IV.4.2. We proceed as in Theorem IV.3.2. We fix (V,TZ, Vs) and we minimize the 
functional v i — > Jls(K,T^-,Vs,v) over the space 

W= |V e (L 2 (0,L;H\lj))) 3 I J v'(S 1 ,S 2 ,s 3 )dS 1 dS 2 = for a.e. s 3 G (0, L)|. 

Through solving simple variational problems (see [14] again), we find that the minimum of this functional is 
obtained for 

dV s Si - S 2 dK dK 



(JV.4.19) 



v(S u S 2 , .)-n 1 = -v\S 1 -±-t + \ 1 — • n 2 + SiS 2 — • ni 
L ds3 I as3 ds3 J 

,c r, s (n d Vs „ „ dTZ Sl-SfdTZ i 

v( Sl ,S 2 , .) • n 2 = -,{s 2 — ■ t - Sl S 2 - ■ n 2 + -^-^ • m} 

/ „ n N ry dVs „ dVs , r, n \ <^ 

U(Si,S 2 ,. -t + 5i— - -m + S 2 —^- n 2 = x (S u S 2 )—-t 
dss dss dss 
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Then the symmetric tensor E (see (IV. 3. 7)) at the minimum is given by 

( -v% 3 (K,V s ) 

, _ tos , i(£ + *)g.t 



(JV.4.20) 



E(^,V 5 ) 



2Va5i S2 ^ -t 



E 33 (^,V S ) / 



where £33(7?., Vs) = -7— ■ t — Si— — • ri2 + S2- — • ni. Upon replacing E by E(K, Vs) in the expression of 
ds 3 as 3 ds 3 

Jls and using (IV. 4. 16), we obtain that the minimum of the functional Jls{U,H, Vs, •) over the space W 
is equal to: 

• if k > 3 



Eh f L (dK \ 2 EI 2 f L (dK 

—Jo w n v + —Jo w ni 



fj,K 
~4~ 



as* 



then, we immediately deduce that the minimum m K of the above quantity is obtained with Uo = IZq = 
and it is given by the minimum of the functional Tls defined by (IV. 4. 18). 

• if k = 3 



Eh 



( 



dK 



2 Jo ^ ds 3 



n 2 



2 EI 2 f L (dK \ 2 [iK f L (dK 

+ —j U- ni ) +- 







f L f\\ 
Jo 



K A t\\j 



+- 



E\cv 



where the relation between Vs, U and Ue is given by (IV. 4. 16). Now, if f(ss) > for a.c. S3 e (0, L), then 
the minimum m 3 of the above quantity is obtained with Uo = Ko = 0. 

We now prove that under the condition (IV. 4. 10) we still have Uo =Ko = 0. To this let (U, K, ,v) (we 
have chosen Vs = 0) be in \5u n and vs be a sequence of admissible deformations given by Theorem (IV. 3.1) 
such that 



J LS {U,K, 0,v) = lim 



S-s-0 



Jjvs) - J (h) 

S 2k 



and 



l|E|| ( L2(j2))3x3 = ||E|| (L 2 (0)) 3x3 = Yua—\\dist(V x vs,SO(3))\\ L 2(p s y 
In view of (IV.4.3), (IV.4.4), (IV.4.6) and (IV.4.8) we obtain 

(f - ^MII/II(lw)) 3 )|I e II(lw- ^ Jls(u,h, o,u). 

Now we choose v as the minimizer of Jls{U,K, , •) over W in the above inequality, it gives 



(4 ~ ^~l w II^H(i 2 (o,i)) 3 
.Eh f L f 



f L ?,( (dK \ 2 (dK \ 2 i K f L (dK \ 

i < 1+2 " 2 >WsH + *(sr-)} + T/ (35-*) 



2 7 \ds 3 



2 Eh 
2 



\ds-- 



■ ni 



4 ./o ^S 3 



/||ftAt|||. 



It follows from the above analysis that the minimum m 3 is obtained for Uo = T^o = 0. In both cases the 
minimum m 3 is given by the minimum of the functional Tls defined by (IV. 4. 18).) □ 
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Remark. Still in the case where k = 3, if / satisfies (IV. 4. 2) without any other assumption, the limit 
deformation is given by Theorem IV. 2. 2. 

We assume k — 3. Let us end this section with an interesting model obtained by superposing general 
forces given by (IV. 1.5) denoted here by (fs,9s) an d the special forces f} described in (IV. 4. 2) satisfying 
the assumptions of this section ((IV. 4. 10) or (IV.4.13)). The corresponding total limit energy in this case is 
then 

\/{U,Tl,Vs,v) e Uun, Jlg{U,TI,Vs,v) = J {^(ME)) 2 + m|||E||| 2 } 

2 L L 

~tlj (/ 9°(S u S 2 ,.)S a dS 1 dS 2 )-(TlAn a )-\Lj\^ f 1 ^ • t. 

Notice that the above energy leads to a linear limit model which couples the inextentional displacement, the 
torsion and the extentional displacement (see also (IV. 4. 24)). 

The analysis of the present subsection and of Subsection IV. 3 permits to state the following theorem. 

Theorem IV.4.3 The functional J LG is the T-limit of ^/^^ in the sense of Theorem IV. 4.1. 

0° 

Let (Uo,lZo,Vs.o,v ) be a minimizer of the functional Jlg over Uh« and define Ue,o € D Ex by 
(IV. 4-16). Then (Uq 7 TZo, Ue,o) is the unique solution of the following minimization problem: 



(/y.4.21) 

where 



m 3 = FLGiUo,Tlo,UEfl) — min Tlg{U,K,Ue), 

(u,n,UE)eVGi in 



VG lin ^{(U,n,U E )e{H 2 (0,L)) 3 x(H 1 (0,L)) 3 xDE X W(0) = R(0) = 0, U E {0)=0, — =KAt}, 
and 



(7V.4.22) 



T lg {U,TI,Ue) 



Eh 



L r 



all 

ds 3 



n 2 



+ - 



EU 



dU t 
ds 3 



, EI 2 


r L 


-dK 


r 


2 J 







-ds 3 


1 


dU 


2- 


2 




2 


ds 3 


2- 







ni 



2 uK 
+ t -r- 



— t] 2 

ds 3 . 



L 2 

-j (Uf + J2J' 9°(SuS 2 ,.)S a det(m |n 2 | ^)dS 1 dS 2 ) U 
~ E £ (/ 9°(Si,S 2 , .)S Q dS!d5 2 ) • (U A n Q ) - | w | £ f 1 ^ • t 



E is the Young's modulus and K is given in Theorem IV. 1.3. Moreover we have 

dU 



(IVA.23) 



dU E ,o 1 
ds 3 ' 2 



ds 3 



2 o _ 

2 + f'' 



and i/ie couple {U ,H ) £ V/ in is £/ie unique solution of the following variational problem: 



(7V.4.24) < 



i>L 2 



ds 3 

2 



rift 

ds 3 



dK 



ds 3 



■ t 



dK 

ds 3 



■ t 



H f L j x &h dU_ 
2 J ds 3 ds 3 



I + 9°(S u S 2 ,.)S a det(n 1 | n 2 | ^)dSidS 2 ) ^ 

+ V / ( / ff°(Si,S 2 , .)S a dSidS 2 ) ■ (n A Da), V(H,ft) G V Kn 
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Indeed the remarks at the end of Sections IV. 2 and IV. 3 are still valid for the above chosen forces. 

V. Solutions of the non-linear minimization problem (IV.2.21) 

The results of this subsection are limited to the case where the curved rod is fixed only on r ,<5 (see 
Subsection II. 2. 4). As a consequence, the other extremity (for S3 = L) is free (or with little change submitted 
to a given load). For these boundary conditions, we replace the minimization problem (IV.2.21) by an integro- 
diffcrcntial equation satisfies by R. To do that, we write the minimization problem (IV.2.21) in terms of the 
unknown R. We denote by G the matrix of (L 2 (0,L)) 3x3 such that 

£ < G,R- I 3 > = £ (\oj\f + J2j9(Si,S 2 ,.)S a det (m | n 2 | ^-)dS l( iS 2 ) ■ (V - M) 

+ E/ (J 9(Si,S2,.)S a dS 1 dS 2 ) -(R-IsK, 

for any (R, V) <G V n lin, where < •, • > is the inner product associated to the Frobenius norm over the space 
M 3 . 

We set 

A3 ={A g (L 2 (0,L)) 3x3 I A T (s 3 ) = -A(s 3 ) for a.c. s 3 e (0,L)} 
US ={R E (i7 1 (0,L)) 3x3 I R(0) = I 3 and for any s 3 G [0,L], R(s 3 ) € 50(3)}. 

Let A be a matrix belonging to A3 and let Ra be the solution of the Cauchy's problem 



(V.1) 



f R A e (i/ 1 (0,i)) 3x3 , 

^-^(s 3 ) = R A (s 3 )A(s 3 ), for a.e. s 3 e (0, L), 

as 3 

[R A (0)=I 3 . 



The map A 1 — > R A is one to one from A3 onto US. An element R e US is associated to the element 

A = R T ^of^ 3 . 

as 3 

Taking into account the definition of G, the minimum to 2 is in fact the minimum of the functional 
(V.2) ^(R) = T y o E/«(^t-Rn a ) +^j o (^n,R„ 2 ) - ^ <G,R I 3 > 

over the closed set US. In terms of A, mi is also the minimum of the functional 

m fl(A, - f £ ± 1. (^t . a A „„)% <f £ (f± „, . r a „ 2 ) 2 - £ < g, a A - u > 

over the space ^3. In view of (V.l), we have 

fL 2 n nj^ r> pL 



iVA) QM = §1 E/ Q (At.„„) 2 + ^^ (Anx-n^-jf 



< G , R A I3 > 
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In what follows we derive the first and the second derivatives of the last term in (V.4). By a standard 
calculation we show that for any matrices A and B in A3 we have 

Ra+b(s 3 ) =Ra(s 3 ) + ( R A (s)B(s)R A (s)ds)R A (s 3 ) 

+ (J 3 J R A (t)B(t)R^(t)R A ( s )B( s )R^( s )d^ s )R A ( S 3)+O(|||B|||3 L2(0ji)) 3 X3 ), 



as the consequence we obtain 



g(A + b) = g(A) + g'(A)(B) + -g"(a)(b, b) + o(\\b\\ 3 {l2{0!L))3x3 ), 



where 



(V.5) 



Q (A) (B) = E [ V I a (At . n a ) (Bt • n Q ) + ^ /" (Am • n 2 ) (Bin ■ n 2 ) 

Jo Q=1 2 Jo 

-jf <G(* 3 ),(jf R A (s)B(s)R A (s)ds)R A (s 3 ) >ds 3 

g"(A)(B,B)=E f V/ Q (Bt-n Q ) 2 + ^ / (Bm • n 2 ) 2 
Jo a=1 Jo 

-2^ <G(« 3 ),(j£ ^ R A (i)B(i)R A (t)R A (s)B(s)R A (s)dtds)R A (s 3 ) >ds 3 . 

In order to explicit the minimum of we simplify the term involving the forces in (A)(B). We have 

<G(a 3 ),(J R A (s)B(s)R A (s)ds)R A (s 3 ) > ds 3 
= jf < G( S3 )R A ( S3 ) , ( jf 3 R A ( s )B( s )R A ( s )d s ) > ds 3 . 

We integrate by parts the right hand side term in the above equality. This gives 

J < G(s 3 ) , ( R A (s)B(s)R A (s)ds)R A (s 3 ) > ds 3 
= L < (/ G(s)R A (s)ds),R A (s 3 )B(s 3 )R A (s 3 ) >ds 3 
= jf <R A (s 3 )(^ G( S )R^(s)ds)R A ( S3 ),B( S3 ) > ds 3 . 

Using the fact that symmetric and antisymmetric matrices are orthogonal for the scalar product < • , • > , we 
finally get for any matrix B e _4 3 



(V.6) 



g'(A)(B)=E^ ^ J Q (At-n a )(Bt-n a ) + ^jf (Am • n 2 ) (Bm • n 2 ) 

-£ < Rl(s 3 )(£ i[G( s )R A ( s ) -R A ( s )G T (s)]d s )R A ( S3 ),B( S3 ) > ds 3 . 



The above derivations allow to prove the following theorem. 
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T dRo 



Theorem V.l. Let (Vo,Ro) be in V n iin and set Ao = R — — . Then (Vo,Ro) is a solution of the 

as 3 

minimization problem (IV. 2. 21) if and only if Ho is a solution of the following integro- differential problem 
A (s3)ni(s 3 ) • n 2 (s 3 ) = 7^ R o [GRq - RoG T ])Ro(s 3 )n 1 ( S3 ) ■ n 2 (s 3 ) 



(V.7) 



Ao(s 3 )t(s 3 ) • ni(« 3 ) = -^R^(s 3 )(^ [GRj - R G T ])R ( S3 )t( S3 ) • m(s 3 ) 
A (s 3 )t(s 3 ) • n 2 (s 3 ) = -2_R^( S3 )(^ [GRj - R G T ])R ( S3 )t( S3 ) • n 2 (s 3 ). 



Moreover, if 

(V-8) ||G|| (i 2 ( o,i))3x3 < -^inf (Eh,EI 2 ,^) 

the solution of the minimization problem (IV. 2. 21) is unique. 

Proof. An element (Vo, Ro) of V„; m is a minimizer of (IV. 2. 21) only if Ao is a minimizer of the functional 
G given by (V.3). Hence, we have £/'(A )(B) = for any B e A 3 . In view of (V.5) and (V.6), the 
antisymmetric matrix Ao satisfies 

E f 7 ° ( A o* ' n «) ( Bt ' n «) + I ( A ° n i ' ( Bn i ' n a) 

Jo Q=1 ' Jo 

= £ < Ro (*a) ( \ [G{s]B% - RGJ] ) R (s 3 ) , B(s 3 ) > ds 3 , VBei 3 - 
This immediately gives (V.7). 

Now we prove that the functional Q admits a unique minimizer, under the assumption (V.8). For any Aei 3 
we get 

lll|2 fn i|2 ii m2 ii m2 1 

ll A ll(L 2 (0,L) 3 >< 3 = 2 \ll At ' ni llL 2 (0,L) + ll At ' n2 lli 2 (0,L) + ll Ani ' n2 lli 2 (0,L) J " 

From the expression (V.5) of Q (A)(B, B) and the above equality we have 
(V.9) 0"(A)(B,B) > \{ inf (EI U EI 2 , ^) - L 3 / 2 \\G\\ {ma ^ 3 }\\B\\l 



2 

L 2 (0,L)) 3x3 - 



As a consequence of the above inequality, if G satisfies (V.8) the functional Q is strictly convex, which insures 
the uniqueness of the minimizer Ao. □ 

Appendix. A few recalls on rotations 

Let V be a matrix belonging to 5*0(3). The matrix V is the matrix of a rotation 7\L a ,e in R 3 where a is 
a unit vector belonging to the axis of the rotation and where 6 belonging to [0,tt] is the angle of rotation 
about this axis. The rotation is written as 

Vx e M 3 , fc a ,e(x) = cos(#)x + (1 - cos(0)) < x, a > a + sin(0) a A x. 

We have 

0\ 2y/2, 



|I 3 -V|||=2V2sin(|)>^. 
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For all t G [0, 1], we denote by W(i) the matrix of the rotation TZ^te- The function t — > W(t) belongs to 
((^([O,!])) 3 * 3 and satisfies 

W(0) = I 3 , W(1) = V, W(t)eSO(3), |||^(*)||| = v / 20<2|||I 3 -V|||, te[0,l]. 



Now, let Uo and Ui be two elements in SO(3). We set 



V = U C T 1 U 1 



and we consider the map 



U(t) = U W(t) for any i G [0, 1], 



where W(t) is defined above. We have built a path U G (C 1 ([0, l])) 3x3 such that 



U(0) = U , U(l) = Ui, U(t)eSO(3), 
dU, 



(*) 



<2|||Ui-U | 



te [0,1]. 



Lemma A. Le£ R 6e in (i/ 1 (0, L)) 3x3 such that R(0) = I3 and swc/i iftai /or any s 3 G [0, L] ifte matrix 
R(s 3 ) belongs to SO(3). There exists a sequence of matrices (Rjv)jveN* satisfying Rat € (VF 1,oo (0, L)) 3x3 , 
Rjv(0) = I3 and for any s 3 G [0,L] i/ie matrix Rat(s 3 ) belongs to 50(3) and moreover 



R 



R strongly in (H l (0,1/)) 



3x3 



dR 

Proof. The matrix A = R T -r — is antisymmetric and belongs to (£ 2 (0, L)) 3x3 . Let (AAr) ngN be a sequence 

US 3 

of antisymmetric matrices such that 



AjvG (C([0,L])) 



3x3 



and 



*-N 



A strongly in (L 2 (0,L)) 



3x3 



Let Rtv (AT G N) be the solution of the Cauchy's problem 

dRjv 



— — — RtvAat 
ds 3 

Rw(0) = I 3 

We have Rjv G (C 1 ([0, L])) 3x3 and for any s 3 G [0,L] the matrix Rjv(s 3 ) belongs to 50(3). From the above 
strong convergence we deduce that 



R 



■N 



R strongly in (id 1 (0,1,)) 



3x3 



□ 
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